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ABSTRACT 

We present Herschel observations of 62 early-type galaxies (ETGs), including 39 galaxies morpho- 
logically classified as SO+SOa and 23 galaxies classified as ellipticals using SPIRE at 250, 350 and 
500 /urn as part of the volume- limited Herschel Reference Survey (HRS). We detect dust emission in 
24% of the ellipticals and 62% of the SOs. The mean temperature of the dust is (T d ) = 23.9 ± 0.8 K, 
warmer than that found for late-type galaxies in the Virgo Cluster. The mean dust mass for the 
entire detected early- type sample is logM^ = 6.1 ± 0.1 M with mean dust-to-stellar mass ratio of 
log ( Aid /M*) = —4.3 ± 0.1. Including the non-detections, these parameters are logM^ = 5.6 ± 0.1 
and log(M,j/M*) = —5.1 ± 0.1 respectively. The average dust-to-stellar mass ratio for the early-type 
sample is fifty times lower, with larger dispersion, than the spiral galaxies observed as part of the 
HRS, and there is an order of magnitude decline in M^/M* between the SOs and ellipticals. We use 
UV and optical photometry to show that virtually all the galaxies lie close to the red sequence yet the 
large number of detections of cool dust, the gas-to-dust ratios and the ratios of far-infrared to radio 
emission all suggest that many ETGs contain a cool interstellar medium similar to that in late-type 
galaxies. We show that the sizes of the dust sources in SOs are much smaller than those in early-type 
spirals and the decrease in the dust-to-stellar mass ratio from early-type spirals to SOs cannot simply 
be explained by an increase in the bulge-to-disk ratio. These results suggest that the disks in SOs 
contain much less dust (and presumably gas) than the disks of early-type spirals and this cannot be 
explained simply by current environmental effects, such as ram-pressure stripping. The wide range in 
the dust-to-stellar mass ratio for ETGs and the lack of a correlation between dust mass and optical 
luminosity suggest that much of the dust in the ETGs detected by Herschel has been acquired as the 
result of interactions, although we show these are unlikely to have had a major effect on the stellar 
masses of the ETGs. The Herschel observations tentatively suggest that in the most massive systems, 
the mass of interstellar medium is unconnected to the evolution of the stellar populations in these 
galaxies. 
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1. INTRODUCTION 

In the standard theoretical view of galactic evolution, 
based on the hierarchical paradigm, galaxies may move 
both ways along the Hubble sequence during their evo- 
lution: ellipticals may form as the result of the merging 
of late-type galaxies (e.g., Cole et al. 2000; De Lucia et 
al. 2006); disks form as the result of the accretion of gas 
on to bulges (Cole et al. 2000). There are also plenty of 
environmental processes that might transform the mor- 
phology of a galaxy, just one example among many being 
ram-pressure stripping in clusters, which would move a 
galaxy toward the early-type sequence by removing the 
interstellar medium (ISM; e.g., Corbelli et al. 2011) and 
thus quenching star formation (Boselli & Gavazzi 2006). 

The observational evidence that galaxy transformation 
is occurring is mixed, especially for galaxies at the cur- 
rent epoch. The fact that the morphology-density rela- 
tion is not just the result of early-type galaxies (ETGs) — 
ellipticals and SOs — being more common in clusters, but 
extends over a wide range of environmental density, is 
evidence that galaxy transformation is not a significant 
process at the current epoch (Dressier 1980). However, 
there is also plenty of observational evidence that merg- 
ing is still important today. In particular, there is ev- 
idence of recent merging or accretion in nearby ellipti- 
cals, primarily from obscuration in the optical (Goud- 
frooij & de Jong 1995), where fossils of mergers in the 
form of dust ripples and dust lanes have been detected 
(Schwcizcr & Seitzer 1992). Furthermore, the discovery 
with the Infrared Astronomical Satellite (IRAS) of the 
ultraluminous infrared galaxies (ULIRGs), which are of- 
ten late-stage mergers with optical profiles similar to the 
profiles of ellipticals (Wright et al. 1990), was persuasive 
evidence that ellipticals are still being formed by merg- 
ers at the current epoch. Nevertheless, it is not clear 
whether all the properties of the eventual merged sys- 
tems would be the same as present-day ellipticals (Naab 
& Ostriker 2009), and many of the distinctive properties 
of ellipticals seem more naturally explained if ellipticals 
form in a relatively short early-period in cosmic history 
(Peebles & Nusser 2010). 

The properties of the ISM in ETGs potentially have 
something to tell us about the evolution of these systems. 
Since late-type galaxies are known to contain much more 
gas and dust than ETGs, any process that moves a galaxy 
along the Hubble sequence has to change simultaneously 
the morphology and the ISM of the galaxy. In this paper, 
we present Herschel observations of continuum emission 
from the dust in a sample of the most massive ETGs in 
the nearby universe, which we show is the most sensitive 
method currently available for estimating the mass of the 
cool ISM in these galaxies. 

Although once thought to be devoid of a cool ISM, op- 
tical absorption studies suggested that more than 50% of 
ellipticals contain some dust (Goudfrooij et al. 1994; van 
Dokkum & Franx 1995; Tran et al. 2001; Ferrarese et al. 
2006) and therefore molecular gas, though the mass of 
dust is uncertain from these works. In the last decade, 
previous far-infrared observations with IRAS, the In- 
frared Space Observatory (ISO), Spitzer, and AKARI 
have shown that dust in ETGs is rather common, with 
dust masses estimated at some 10-100 times larger than 
those derived from the optical extinction measurements 



(e.g., Bregman et al. 1998; Ferrari et al. 2002; Kcnnicutt 
et al. 2003; Pahre et al. 2004; Xilouris et al. 2004; Temi 
et al. 2004, 2007a; Panuzzo et al. 2007; Kaneda et al. 
2008; Young et al. 2009). Although these observations 
have probed the peak of the spectral energy distributions 
(SEDs) of the brightest giant ellipticals and SOs, they are 
hampered by poor resolution and/or lack of long wave- 
length coverage. 

Our current ideas about the life cycle of dust in el- 
lipticals suggest that it may be possible for us to test 
the origin of the dust, particularly whether it is provided 
by internal (via stellar mass loss) or external processes 
(e.g., fueled by mergers). As long as dust is not intro- 
duced from elsewhere, we essentially have a stellar sys- 
tem in equilibrium: dust is produced in the atmospheres 
of evolved stars (Athey et al. 2002) and destroyed via 
sputtering in the hot gas (indeed Bressan et al. 2006 
detect silicate features in elliptical galaxies known to be 
produced in circumstellar envelopes). The dust life cy- 
cle in ellipticals may therefore be simpler than that in 
late-types for the internal origin scenario (Tsai & Math- 
ews 1996), and the mass of dust predicted to exist in 
the steady state through production in evolved stars and 
destruction through sputtering is approximately 10 5 M Q 
(Forbes 1991; Goudfrooij & de Jong 1995). 

An alternative origin for dust in ETGs is accreted ma- 
terial either due to a merger or from a tidally-interacting 
companion: Centaurus A has long been known to har- 
bor a dusty disk thought to be formed by a merger (re- 
cently revealed in the submillimetcr (submm) — Leeuw 
et al. 2002; Auld et al. 2011; T. Parkin et al. 2011, sub- 
mitted). Forbes (1991) and Temi et al. (2004, 2007a) 
found no correlation between the dust emission and op- 
tical starlight using ISO observations of massive ETGs, 
suggesting that the dust has an external rather than an 
internal origin. In 2010, Gomez et al. used Herschel ob- 
servations of the elliptical galaxy M86 to reveal 10 6 M© 
of cold dust coincident with material stripped from the 
nearby spiral NGC4438 (Cortese et al. 2010a). A similar 
process has been seen with the "displaced ISM" of the 
elliptical NGC3077 (Walter et al. 2011). 

The unprecedented resolution and sensitivity of the re- 
cently launched Herschel Space Observatory (Pilbratt et 
al. 2010), combined with the wavelength coverage of the 
instruments PACS (Poglitsch et al. 2010) and SPIRE 
(Griffin et al. 2010) from 70 to 500 /im, allow us to ad- 
dress long-standing issues such as the origin and quantity 
of dust in ETGs. Herschel provides us with an unbiased 
view of the interstellar dust (and therefore the interstel- 
lar medium) ; we are now sensitive to the total dust mass 
in galaxies rather than only to those galaxies which had 
enough warm dust to be detected by IRAS and Spitzer 
(e.g., Dunne et al. 2011). Here, we use observations from 
the guaranteed time project, the Herschel Reference Sur- 
vey (HRS; Boselli et al. 2010b) of 62 ETGs including 23 
ellipticals. The HRS is a study of dust in nearby galax- 
ies, crucial for calibrating blind surveys of galaxies with 
Herschel at high redshifts, e.g., HerMES (Oliver et al. 
2010) and the Herschel- ATLAS (H- ATLAS; Eales et al. 
2010). 

Cortese et al. (2011b, hereafter C12) recently used the 
HRS sample of ~300 nearby galaxies to obtain dust scal- 
ing relations, finding that the dust is tightly coupled to 
the atomic mass and therefore to the cold ISM. C12 mea- 
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sured the dust-to-stellar-mass ratio of galaxies with dif- 
ferent morphological types and environment, finding that 
the ETGs contain less dust mass per unit of stellar mass 
than the late-type galaxies. In this work, we further ex- 
plore the dust properties of the ETGs in the HRS and 
measure more precisely the variation of dust content with 
Hubble type and the origin of this variation. In Section [2] 
we introduce the sample and describe the data reduction 
techniques. The results are given in Section [3] and are 
used to discuss our understanding of the dust content of 
early- and late-types in Section l4~T1 and the link between 
dust and the hot interstellar medium in Section l4~2l We 
discuss the implications for the evolutionary history of 
ETGs in Section [5l The conclusions are presented in 
Section [6] 

2. THE SAMPLE 

The HRS is a volume-limited sample (15Mpc < d < 
25 Mpc) with Herschel SPIRE of 322 galaxies selected by 
the if- band magnitude (a proxy for stellar mass). The 
ETGs in the HRS have K < 8.7, and so our sample con- 
sists of the ETGs in this volume of space with the high- 
est stellar masses; contamination from Galactic cirrus is 
minimized due to the high galactic latitudes of the sam- 
ple (see Boselli et al. 2010b for full details). The sample 
includes different morphological types with 260 late-type 
galaxies and, in the original catalog listed in Boselli et 
al. (2010b), 64 early types. However, after closer in- 
spection of the optical images, we reclassified the galax- 
ies NGC4438 (HRS163), NGC4457 (HRS173), NGC4691 
(HRS256), and NGC5701 (HRS322) as late types (from 
SO/a to Sb). We changed two classifications in the other 
direction, with the Sb galaxies NGC4179 (HRS90) and 
IC3510 (HRS202) reclassified to SO and dwarf elliptical, 
respectively. The revised catalog therefore has 62 ETGs, 
39 of which are designated SO+SOa galaxies and 23 are 
ellipticals (Table Q} though five of the latter set are mor- 
phologically classified as E/S0. The morphological clas- 
sifications are taken from the Virgo Cluster Catalogue 
(Binggeli et al. 1985) or are our own if another was not 
available (Boselli et al. 2010b). The HRS contains galax- 
ies from a range of environments, including isolated field 
galaxies, pairs and galaxies in the heart of the Virgo Clus- 
ter. 

The HRS observations were carried out with SPIRE 
in scan-map mode with a scan speed of 30 arcsec s . 
We chose the size of our maps so that they would be at 
least as large as the optical disk of each galaxy, defined 
as the area of the galaxy within an optical isophote with 
a _B-band brightness of 25 mag arcsec -2 (D(25); Boselli 
et al. 2010b). In practice, our maps were either 4x4, 
8x8, or 12x12 arcmin 2 , and in most cases are larger than 
13(25). For each ETG, we made eight pairs of orthogonal 
scans, resulting in a la instrumental noise in each pixel of 
0.35, 0.20 and 0.11 MJy/sr for the 250, 350, and 500 ^m, 
respectively (integration time of 1199s, 3102s, and 4948s 
for each map size). We deliberately observed the ETGs 
in the HRS for at least twice as long as for the late-type 
galaxies (three pairs of orthogonal scans) because the 
former are known to contain much less dust on average. 

There are a number of sources in the HRS which 
overlap with galaxies observed as part of the open-time 
key project the Herschel Virgo Cluster Survey (HeViCS; 
Davies et al. 2010, 2012), thus the two surveys have a 



data-sharing agreement for these galaxies. Here we also 
include data of the duplicate 19 ETGs formally observed 
as part of HeViCS with both PACS and SPIRE. The 
SPIRE observations for HeViCS were made in parallel- 
scan-map mode with a scan speed of 60 arcsec s . To 
cover the Virgo Cluster, four fields with size 4° x 4° were 
observed with eight cross-scans (see Davies et al. 2010, 
2012 for a complete description). Note that the dust 
content of all of the ETGs in the Virgo Cluster from 
HeViCS will be presented in a complementary paper (S. 
di Serego Alighieri et al., in preparation). 

2.1. Data Reduction and Flux Extraction 

The HRS and HeViCS SPIRE data were reduced us- 
ing similar pipeline procedures. The SPIRE data were 
processed up to Level- 1 with a custom script adapted 
from the official pipeline (P 'OF '5 -pipeline. py, dated 8 
Jun 2010) as provided by the SPIRE Instrument Con- 
trol Centre (ICC^B Our custom Jython script was 
run in the Herschel Interactive Processing Environment 
(HIPE; Ott 2010) with the continuous integration build 
number: 4.0.1367. For both surveys we use an op- 
timized deglitcher setting instead of applying the ICC 
default settings. For the HeViCS data the sigmaKap- 
paDeglitcher module was used, while for the HRS we 
applied the waveletDeglitcher; this module was ad- 
justed to mask the sample following a glitch. For the 
HRS, after the flux calibration was applied, an additional 
pass with the waveletDeglitcher was run, as this was 
found to significantly improve the glitch-removal process. 
Furthermore, we did not run the pipeline default temper- 
ature drift correction or the median baseline subtraction. 
Instead we use a custom method (BriGAdE; M. W. L. 
Smith et al., in preparation; see also L. Ciesla et al., in 
preparation) to remove the temperature drift and bring 
all bolometers to the same level. 

Our final SPIRE maps were created using the naive 
mapper provided in the standard pipeline with pixel sizes 
of 6", 8", and 12" at 250, 350, and 500 ^m respectively. 
The FWHM of the SPIRE beams for this pixel scale are 
18.2", 24.5", and 36.0" at 250, 350, and 500 ^m, respec- 
tively (Swinyard et al. 2010). In addition, the 350 /mi 
data are multiplied by 1.0067 to update our flux densi- 
ties to the latest v7 calibration product. The calibra- 
tion uncertainty is a combination of the 5% error due 
to correlated errors between bands and the 2% random 
uncertainty; these values are added linearly instead of in 
quadrature^. 

The Herschel PACS 100 and 160 /mr data taken as part 
of HeViCS were reduced using the standard pipeline (see 
Davies et al. 2012). In brief, we used a two step deglitch- 
ing process, first with the standard deglitcher and second 
using one based on sigma clipping. We masked bright 
sources and a high-pass filter was used to reduce 1// 
noise. The orthogonal scans were combined and the naive 
mapper was used to create the final maps. Unlike Davies 
et al. (2012) we use images created using the full HeViCS 
data set with eight parallel-mode scans per tile. The 

18 See 'The SPIRE Analogue Signal Chain and Photometer De- 
tector Data Processing Pipeline' (Griffin et al. 2008 or Dowell et 
al. 2010) for a more detailed description of the pipeline and a list 
of the individual modules. 

19 SPIRE Observer's Manual (2011) 
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The Sample 



HRS 


Other Name 




R.A. 




Decl. 


Type 


D 


D(25) 
(') 


L B 


L K 


Membership 








(h 


m 


s) 


(° ' ") 




(Mpc) 


(L@) 


(L Q ) 




3 


NGC 


3226 


10 


23 


27 


,01 r 


+19 53 54. 7 r 


E2:pec;LINER;Sy3 


16.7 


3.16 


10.12 


10.59 


Leo CI. 


7 


NGC 


3245 


10 


27 


18. 


40 r 


+28 30 26. 3 r 


SA(r)0:?;Hll; LINER 


18.8 


3.24 


10.08 


10.47 


Leo CI. 


14 


NGC 


3301 


10 


36 


56 


.04 


+21 52 55.7 


/ t^i i \ o T^t / \ n / 

(R )SB(rs)0/a 


19.2 


3.55 


9.84 


10.47 


Leo CI. 


22 


NGC 


3414 


10 


51 


16. 


19 r 


+27 58 30. 2 r 


On T TAT T - ^ T"i 

SO pcc;LlNER 


20.2 


3.55 


9.98 


10.73 


Leo CI. 


43 


NGC 


3608 


11 


16 


58 


.96 


+18 08 54.9 


E2;LINER: 


15.8 


3.16 


10.11 


10.70 


Leo CI. 


45 


NGC 


3619 


11 


19 


21 


51 r 


+57 45 28. 3 r 


/ T^i \ o A / \ n I 

(R)SA(s)0+: 


22.1 


2.69 


9.87 


10.57 


Ursa Major CI. 


46 


NGC 


3626 


11 


20 


03 


80 r 


+18 21 24. 3 r 


/ T~i \ O A / \ n i 

(R)SA(rs)0+ 


21.3 


2.69 


10.12 


10.70 


Leo CI. 


49 


NGC 


3640 


11 


21 


06. 


85 


+03 14 05.4 


E3 


17.9 


3.98 


10.43 


11.06 


Leo CI. 


71 


NGC 


3945 


11 


53 


13 


61 r 


+60 40 32. 3 r 


O T~» / \ n i T TTVTT IT* 

SB(rs)0+;LINER 


18.0 


5.25 


10.01 


11.06 


Ursa Major CI. 


87 


NGC 


4124 


12 


08 


09 


.64 


+10 22 43.4 


SA(r)0+ 
S0 a 


17.0 


4.10 


9.70 


10.38 


Virgo Out. 


90 


NGC 


4179 


12 


12 


52 


.11 


+01 17 58.9 


17.0 


3.80 


9.89 


10.60 


Virgo Out. 


93 


NGC 


4203 


12 


15 


05 


06 r 


+33 11 50. 2 r 


SAB0-:;LINER;Sy3 


15.6 


3.39 


9.89 


10.73 


Coma I CI. 


101 


NGC 


4251 


"1 o 

Iz 


18 


08 


.31 


+28 10 31.1 


SB0? sp 


15.3 


3.63 


13.86 


10.54 


Coma I CI. 


105 


NGC 


4262 


1 o 
Iz 




30 


58 


1 "1 A P" O O C\ o 

+14 52 39.8 


&B(s)0-! 


17.0 


1.87 


9.69 


10.43 


Virgo A 


123 


NGC 


4324 


1 o 
Iz 


nq 
ZO 


06 


.18 


+05 15 01.5 


SA(r)0+ 


17.0 


3.52 


9.65 


10.38 


Virgo S CI. 


125 


NGC 


4339 


"1 o 
Iz 


Zo 


34. 


,94 


+06 04 54.2 


E0;Sy2 


23.0 


2.31 


9.71 


10.34 


Virgo B 


126 


NGC 


4340 


12 


23 


35 


31 


i 1 c ao 1 n o 

+ 16 43 19.9 


bB(r)0+ 


17.0 


3.60 


9.81 


10.44 


Virgo A 


129 


NGC 


4350 


"1 o 
Iz 


Zo 


57. 


81 


+ 16 41 36.1 


OAn Al l: 

bAU;Abs. line 


17.0 


3.20 


9.91 


10.64 


Virgo A 


135 


NGC 


4365 


12 


24 


28 


23 


i ii7 1 n no -1 

+07 19 03.1 


E3 


23.0 


8.73 


10.34 


11.26 


Virgo B 


137 


NGC 


4371 


1 o 
Iz 


OA 
Z 1 


55. 


.43 


1 "1 "1 A n 1 P" A 

+ 11 42 15.4 


SB(r)0+ 


17.0 


5.10 


9.92 


10.68 


Virgo A 


138 


NGC 4374, M84 


12 


25 


03 


.78 


i 1 n ct i 1 i 1 

+ 12 53 13.1 


tt't t T7 | Tur~' c n 

hjl:Lii;K,G;oyz 


17.0 


10.07 


10.57 


11.26 


Virgo A 


150 


NGC 4406, M86 


Iz 


zo 


11. 


74 


i 1 n r* ac* A 

+ 12 56 46.4 


bv(S)/hj6 


17.0 


11.37 


10.66 


11.31 


Virgo A 


155 


NGC 4417 


12 


26 


50 


.62 


i r\c\ o r no n 

+09 35 03.0 


SB0: s 


23.0 


3.60 


10.08 


10.77 


Virgo B 


161 


NGC 4429 


12 


27 


26. 


56 


+11 06 27.1 


OA/ \n i T TMTTin TTt, 

bA(r)0+;LlNER;TIII 


17.0 


8.12 


10.23 


11.06 


Virgo A 


162 


NGC 4435 


12 


27 


40. 


50 r 


i "in n a a A t*y 

+ 13 04 44.5 


SB(s)0;LINER;Hn 


17.0 


2.92 


10.05 


10.83 


Virgo A 


166 


NGC 4442 


12 


28 


03 


.89 


+09 48 13.0 


SB(s)0 


23.0 


5.05 


10.36 


11.12 


Virgo B 


174 


NGC 4459 


12 


29 


00 


04 r 


+13 58 42. 2 r 


OA / \ n I TT.. T T "\T T — 1 T"> 

bA(r)0+;Hll; LINER 


17.0 


3.36 


10.07 


10.91 


Virgo A 


175 


NGC 4461 


12 


29 


03 


.01 


+13 11 01.5 


SB(s)0+: 


17.0 


3.52 


9.90 


10.57 


Virgo A 


176 


NGC 4469 


12 


29 


28 


.03 


+08 44 59.7 


SB(s)0/a? sp 


23.0 


4.33 


9.97 


10.82 


Virgo B 


178 


NGC 4472, M49 


12 


29 


46. 


76 


+08 00 01.7 


E2/S0;Sy2 


17.0 


10.25 


10.90 


11.59 


Virgo S CI. 


179 


NGC 4473 


12 


29 


48 


.87 


+13 25 45.7 


E5 


17.0 


4.04 


10.15 


10.90 


Virgo A 


180 


NGC 4477 


12 


30 


02 


.17 


+13 38 11.2 


SB(s)0:?;Sy2 


17.0 


3.60 


10.12 


13.77 


Virgo A 


181 


NGC 4478 


12 


30 


17 


.42 


+12 19 42.8 


E2 


17.0 


1.89 


9.89 


10.41 


Virgo A 


183 


NGC 4486, M87 


12 


30 


49. 


.42 


+12 23 28.0 


E+0-1 pcc;NLRG;Sy 


17.0 


11.00 


10.85 


11.43 


Virgo A 


186 


NGC 4494 


12 


31 


24. 


03 


i o et a c on n 

+25 46 29.9 


El-2;Sy 


18.7 


4.79 


10.62 


11.20 


Coma 1 CI. 


200 


NGC 4526 


12 


34 


03 


03 r 


+07 41 57. 3 r 


O A T^i / \ n 

SAB(s)0: 
dE a 


17.0 


7.00 


10.41 


11.18 


Virgo S CI. 


202 


IC 3510 


12 


34 


19. 


33 


+11 04 17.7 


17.0 


1.10 


8.70 


9.20 


Virgo A 


209 


NGC 4546 


12 


35 


29 


51 


-03 47 35.5 


SB(s)0-: 


15.0 


3.31 


10.05 


10.82 


Virgo Out. 


210 


NGC 4550 


12 


35 


30 


61 


i -i n i o 1 p" a 

+12 13 15.4 


OT>n O T TATT1T) 

bBO: sy,L!NER 


17.0 


3.95 


9.66 


10.30 


Virgo A 


211 


NGC 4552, M89 


12 


35 


39 


.88 


i "in on n "i t 

+12 33 21.7 


T7t T TMTTin TTTT O »~l 

E;LlNER;Hll;by2 


17.0 


7.23 


10.29 


11.06 


Virgo A 


214 


NGC 4564 


12 


36 


26. 


99 


i "i "i »~i Ai r* 

+11 26 21.5 


E6 


17.0 


4.33 


9.86 


10.58 


Virgo A 


218 


NGC 4570 


12 


36 


53 


.40 


i r\ry -i a a o n 

+07 14 48.0 


on / ^\ /T7"v 

bO(7)/E7 


17.0 


3.52 


9.96 


10.70 


Virgo S CI. 


219 


NGC 4578 


iz 


O ( 


30 


55 


inn o o "i o A 

+09 33 18.4 


SA(r)0: 


17.0 


3.77 


9.71 


10.41 


Virgo E CI. 


231 


NGC 4596 


12 


39 


55. 


.94 


+10 10 33.9 


bB(r)0+;LINER: 


17.0 


4.76 


10.08 


10.79 


Virgo E CI. 


234 


NGC 4608 


12 


A 1 

41 


13 


.29 


+10 09 20.9 


bB(r)0 


17.0 


4.30 


9.83 


10.51 


Virgo E CI. 


235 


NGC 4612 


12 


A 1 


32 


.76 


i nT "i o p" o n 

+07 18 53.2 


(R)bAB0 


17.0 


2.16 


9.82 


10.35 


Virgo S CI. 


236 


NGC 4621, M59 


12 


42 


02 


.32 


ill o o a o n 

+11 38 48.9 


E5 


17.0 


7.67 


10.32 


11.05 


Virgo E CI. 


240 


NGC 4638 


12 


42 


47. 


43 


i "i "i n/" 1 o n n 

+11 26 32.9 


bO- 


17.0 


2.01 


9.82 


10.49 


Virgo E CI. 


241 


NGC 4636 


12 


42 


49 


.87 


i nn a "i "i r* n 

+02 41 16.0 


T7t /On /"I T TTVTTT'Tl O O 

E/bO/l;LlNER;by3 


17.0 


9.63 


10.51 


11.18 


Virgo S CI. 


243 


NGC 4643 


12 


43 


20 


14 


i n "i p" o a o 1 

+01 58 42.1 


OT> / \r\ 1 T TTvTTT^ Ti /TTt, 

bB(rs)0/a;LlNER/Hn 


17.0 


3.00 


9.98 


10.81 


Virgo Out. 


245 


NGC 4649, M60 


12 


43 


40 


.01 


+11 33 09.4 


E2 


17.0 


5.10 


10.73 


11.46 


Virgo E CI. 


248 


NGC 


4660 


12 


44 


31 


.97 


+11 11 25.9 


E5 


17.0 


1.89 


9.70 


10.47 


Virgo E CI. 


250 


NGC 


4665 


12 


45 


05 


.96 


+03 03 20.5 


bB(s)0/a 


17.0 


4.50 


10.04 


10.80 


Virgo Out. 


253 


NGC 


4684 


12 


47 


17. 


.52 


-02 43 38.6 


bB(r)0+;Hll 


21.3 


2.88 


9.82 


10.42 


Virgo Out. 


258 


NGC 


4697 


12 


48 


35 


.91 


-05 48 03.1 


E6;AGN 


17.7 


7.24 


10.55 


11.16 


Virgo Out. 


260 


NGC 


4710 


12 


49 


38. 


93 r 


+15 09 59.1 r 


bA(r)0+? sp:Hn 


17.0 


4.30 


9.91 


10.74 


Virgo Out. 


269 


NGC 


4754 


12 


52 


17. 


.56 


+11 18 49.2 


SB(r)0-: 


17.0 


5.03 


10.05 


10.81 


Virgo E CI. 


272 


NGC 


4762 


12 


52 


56 


.05 


+11 13 50.9 


bB(r)0 sp;LINER 


17.0 


8.70 


10.21 


10.85 


Virgo E CI. 


286 


NGC 


4866 


12 


59 


27. 


14 


+14 10 15.8 


bA(r)0+;LINER 


17.0 


6.00 


9.84 


10.60 


Virgo Out. 


296 


NGC 


5273 


13 


42 


08. 


36 r 


+35 39 15. l r 


SA(s)0;Syl.5 


15.2 


2.75 


9.54 


10.21 


Canes Vcn. bpur 


312 


NGC 


5576 


14 


21 


03 


.68 


+03 16 15.6 


E3 


21.2 


3.55 


10.16 


10.89 


Virgo-Libra CI. 


316 


NGC 


5638 


14 


29 


40. 


.39 


+03 14 00.2 


El 


23.9 


2.69 


10.09 


10.72 


Virgo-Libra CI. 



Notes. — The columns arc as follows. Column 1: the number of the source in the Herschel Reference Survey (Bosclli ct al. 2010b). Column 
2: other common names for the galaxy. Columns 3 and 4: the right ascension and declination of the galaxy (J2000). This is taken from 
the NASA Extragalactic Database (NED), except positions with a superscript r arc radio positions taken from Wrobcl (1991) and Filho et 
al. (2006). Column 5: morphological type taken from Virgo Cluster Catalogue (Binggcli ct al. 1985) or from our own classification. The 
AGN/Seyfert/LINER classifications were taken from Ho ct al. (1997), Schmitt (2001), and Vcron-Cctty & Veron (2006). Column 6: distance 
in Mpc. We assume all objects in Virgo, including those in the outskirts of Virgo, arc at a distance of 17 Mpc, except for those in Cloud B, 
which wc assume are at a distance of 23 Mpc. For the other galaxies we have calculated distances from the heliocentric velocity using a Hubble 
constant of 70km s _1 Mpc -1 . Column 7: optical isophotal diameter (25 magarcscc -2 ). Column 8: total i^-band luminosity measured using 
Ks tot from 2MASS (Skrutskic et al. 2006). Column 9: cluster or cloud membership from Gavazzi & Bosclli (1999) for Virgo, otherwise from 
Tully (1988) or Nolthcnius (1993) wherever available, or failing that our own estimate. 
a These sources were morphologically classified by hand using optical images. 
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FWHM beam sizes are approximately 9" and 13" with 
pixel sizes of 3.2" and 6.4" for the 100 and 160 fim bands, 
respectiveljH Davies et al. (2012) measured differences 
in global flux densities of up to 20% between sets of cross- 
scans, we therefore choose this as a conservative error for 
our flux estimates as it dominates over other sources of 
uncertainty, e.g., the calibratior£3- 

The 250 /.tm Herschel SPIRE maps for the detected 
S0+S0a and all of the elliptical galaxies in our sample 
are shown in Figures [1] and [5J Most of the objects vis- 
ible on the images are background galaxies, creating a 
potential problem for determining whether our targets 
are detected. Fortunately, these background galaxies are 
almost always unresolved, whereas virtually all our de- 
tections are extended. The full description of the flux 
extraction process and complete photometry information 
for all the HRS galaxies including upper limits on the 
non-detections are provided in L. Ciesla et al. (in prepa- 
ration). The upper limits for the non-detections are es- 
timated using circular apertures with size 0.3 x D25 for 
ellipticals and 0.8 x D25 for SOs (chosen as a conservative 
limit based on the extent of 250 /xm emission seen in the 
detected sample). 

We tested the effect of using different definitions of 
the upper limit on the analysis in this paper. First, 
we used simply the statistical noise on the map, tak- 
ing into account the instrumental noise and confusion 
to give <7 rms . Second, we included the additional uncer- 
tainty in the photometry measurement due to the back- 
ground level (cr s k y ; L. Ciesla et al., in preparation), which 
arises mainly from cirrus contamination. For this work, 
we use the photometry upper limit as in L. Ciesla et 
al. (in preparation), with a — 3(cr rms + Usky) which is 
the most conservative limit on the 250 fim flux for our 
sources. Note that all of the results in this paper are 
also valid if we use 3er rms or 5cr rms (as a typical 3cr or 5cr 
detection limit). 

For the HeViCS galaxies in our sample the PACS 
flux densities we used aperture photometry with suitable 
background regions selected around the source (the as- 
sumed uncertainty is outlined in Scction [2~T|) . We present 
the PACS fluxes in Tabled] 

2.2. Data at other wavelengths 

In addition to the Herschel data, we used 70-160 /im 
data from Spitzer (Kennicutt et al. 2003), reprocessed us- 
ing the techniques described in Bendo et al. (2010b) and 
presented in Bendo et al. (2011, submitted). The cali- 
bration uncertainties were assumed to be 5% at 70 /im 
and 12% at 160 /im (Gordon et al. 2007; Stansberry et 
al. 2007). We use IRAS 60 and 100 /im measurements 
originally presented in Knapp et al. (1989) but modi- 
fied in a private communication to NASA Extragalactic 
Database (NED, Table [2]). The calibration uncertainties 
for IRAS were assumed to be 13% at 60 /*m and 16% at 
100 /im (e.g., Verter & Rickard 1998). Four of our ellip- 
tical galaxies were also detected by ISO and we use the 
fluxes given in Temi et al. (2004). 

Optical photometry for the sample was obtained from 
the Sloan Digital Sky Survey (SDSS) DR7 (Abazajian 
et al. 2009) database. We estimate stellar masses for 

20 PACS Photometer Point Spread Function Document (2010) 

21 PACS Photometer Calibration Document (2011) 



the sample from the z-band luminosities using the re- 
lationship between stellar mass and galaxy color for a 
Chabrier initial mass function (Zibetti et al. 2009). The 
full description of this method applied to the entire HRS 
sample can be found in Cortese et al. (2011) and C12. 
NUV photometry is available from the Galaxy Evolution 
Explorer {GALEX; Martin et al. 2005) from the GR6 
data release (see Boselli et al. 2011; Cortese et al. 2011). 

X-ray luminosities for 38 of our sources were obtained 
from the catalogs of O'Sullivan et al. (2001) and Pelle- 
grini (2010). A further 57 of our sources are included 
in the ATLAS 3D sample of ETGs (Cappellari et al. 
2011a). ATLAS 3D includes molecular hydrogen masses 
estimated from the CO(1-0) intensity (I(CO); Young et 
al. 2011) with a conversion factor from H2 to CO of 
AT(H 2 )//(CO) = 3 x 10 20 cm 2 (Kkms^ 1 )- 1 (e.g., Young 
& Scoville 1982). The CO data only exist for the cen- 
tral region of the galaxies (i.e., within 30") and would be 
an underestimate of the total molecular gas if the real 
gas distribution is extended, though Young et al. find 
that there is no strong evidence for extended emission 
in any of the galaxies except for NGC4649. Atomic Hi 
masses and upper limits were obtained for 48 sources 
from the GOLDMINE database (Gavazzi et al. 2003), 
the ALFALFA survey (Haynes et al. 2011) or Noorder- 
meer et al. (2005) and Springob et al. (2005). The CO 
and Hi masses are presented in Table |U In all cases, 
these have been corrected to the distances assumed for 
the ETGs in this paper. 

3. RESULTS 

In Figure [1] we show the detected SO+SOa galaxies 
and all of the ellipticals in Figure [5] the sources show 
a variety of morphologies, from small blobs to spectac- 
ular ring structures. The difference in the appearance 
of the S0+S0a and ellipticals at 250 /im is very clear in 
these figures. The detected SO+SOa sources are brighter 
than their elliptical counterparts and the FIR emission 
is more extended within the D(25) isophotes. In general, 
this appears to be a trend as we move along the Hubble 
sequence: an aperture of size 0.8 x D(25) is sufficient 
to include the submm emission from the ETGs, yet L. 
Ciesla et al. (in preparation) require an aperture size 
of 1.4 x £)(25) to incorporate all the emission from the 
late-type galaxies. This effect was also noticed by Bendo 
et al. (2007) and Munoz-Mateos et al. (2009). 

The difference between late-types and early-types 
when viewed in the submm and optical wavebands is 
quite startling. For example, at optical wavelengths, the 
giant elliptical NGC4649 (HRS245) is extremely bright 
compared to the smaller, fainter spiral galaxy to the 
northwest (NGC4674). At 250 ^m (Figure© however, 
the situation is reversed, with the spiral galaxy extremely 
bright in the FIR, but no sign of the elliptical. In Ap- 
pendix [A] we provide brief notes on the elliptical galaxies 
that are detected by Herschel. 

3.1. Detection Rates 

We detect 31 out of the 62 ETGs, including 24 
SO+SOa galaxies (62%) and 7 ellipticals (30%, Tabled. 
We detect all three of the ETGs classified as 'pecu- 
liar' (two ellipticals, M87 and NGC3226, and one SO, 
NGC3414). Our detection rates are higher than pub- 
lished in an IRAS study of ETGs (Bregman et al. 
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Figure 1. SDSS gri three-color maps and 250 fim Herschel SPIRE maps of the 24 detected SO+SOa galaxies from the Herschel Reference 
Survey. North is up and east is left. The white aperture shows the elliptical optical B-band isophotal diameter (D(25) at 25 mag arcsec - 2 — 
Table [TJ. The region shown is 8.4' X 8.4'. Red crosses show the location of the VLA FIRST radio detection. The 250 (im beam is shown 
in the corner of (a) NGC3301 and (b) NGC4469. 
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Figure 1. (Continued) 
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Figure 2. SDSS gri three-color maps and 250- fim Herschel SPIRE maps of the entire sample of 23 elliptical galaxies from the HRS and 
HeViCS programs (sources denoted with a (*) indicate that it has been detected by Herschel). North is up and east is left. The white 
aperture shows the optical £?-band isophotal diameter (D(25) at 25 mag arcsec -2 — Table[T}- The boxes are 9.6' X 9.6'. Red crosses show 
the location of the VLA FIRST radio detection. The 250 /xm beam is shown in the corner of (a) NGC3226 and (b) NGC4552. 
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Table 2 

Herschel PACS and IRAS Fluxes for Detected Galaxies. 



HRS Name 


F100 


F160 


F100 


F60 




fJv) 


(Jvl 


(Jvl 


\°y j 


7 






3.97±0.64 


2.09±0.28 


14 






0.92±0.16 


0.48±0.07 


22 






0.56±0.21 


0.25±0.04 


45 






1.83±0.37 


0.38±0.06 


71 






1.36±0.24 


0.26±0.04 


87 






1.75±0.29 


0.43±0.07 


93 






2.16±0.36 


0.59±0.08 


101 








0.12±0.04 


105 




2.71±0.54 


0.39±0.14 


0.18±0.04 


123 


1.80±0.6 


2.42±0.48 


1.99±0.32 


0.41±0.07 


126 






0.37±0.09 


0.09±0.03 


129 






1.09±0.19 


0.36±0.06 


138 


0.78±0.16 


0.46±0.09 


1.16±0.22 


0.50±0.07 


150 






0.33±0.09 


0.11±0.04 


161 


4.74±0.95 


4.54±0.91 


5.15±0.83 


1.56±0.21 


162 


4.48±0.90 


3.75±0.75 






174 


4.59±0.92 


3.78±0.75 


4.82±0.78 


1.87±0.25 


176 


3.33±0.67 


3.89±0.78 


3.40±0.57 


1.02±0.14 


180 


0.94±0.19 


0.97±0.19 


1.41±0.25 


0.57±0.09 


183 


0.38±0.08 


0.60±0.12 


0.41±0.12 


0.39±0.07 


186 








0.19±0.06 


200 


15.81±3.16 


15.53±3.11 


17.10±2.74 


5.56±0.72 


209 






0.89±0.2o 


0.26±0.06 


210 


0.28±0.06 


0.27±0.05 


0.25±0.10 


0.14±0.04 


211 






0.53±0.10 


0.16±0.05 


231 






0.75±0.13 


0.40±0.06 


243 






2.06±0.34 


0.62±0.09 


253 






2.15±0.36 


1.27±0.17 


258 






1.24±0.21 


0.46±0.06 


260 






14.79±2.37 


5.73±0.75 


286 






1.02±0.23 


0.15±0.06 


296 






1.56±0.28 


0.90±0.12 


312 






0.21±0.28 


0.09±0.03 


316 






0.45±0.14 




Notes. — 


The Spitzer, and Herschel photometry is presented in Temi 



ct al. (2004), G. Bcndo ct al. (2011, submitted), and L. Cicsla ct al. 
(in preparation), respectively. The columns are as follows. Column 1: 
the number of the source in the Herschel Reference Survey (Bosclli ct 
al. 2010b). Columns 2 and 3: Herschel PACS fluxes at 100 and 160 fim. 
Columns 4 and 5: IRAS fluxes at 100 and 60 fj,m, 

1998), which detected 12%-17%. Although the IRAS 
study excluded peculiar galaxies and active galactic nu- 
clei (AGNs; with approximately half their sample at dis- 
tances larger than the HRS), the small number of these 
sources in our sample suggests that the most likely ex- 
planation of the large difference in detection rates is that 
Herschel is more sensitive to cold dust than IRAS. This 
is supported by the higher mean dust temperature re- 
ported by Bregman et al. (1998). Temi et al. (2004) de- 
tected a higher fraction of ETGs (41% of ellipticals and 
79% of SOs) with ISO, but their sample contains pecu- 
liar galaxies and giant ellipticals. The detection rates for 
the SPIRE and multiwavelength data for our sample are 
listed in Table [3] these rates are significanthP 2 ! different 
for SOs and ellipticals at the 96% level. Our overall detec- 
tion rate of 50% for the ETGs is higher than the 22% ob- 
tained by the ATLAS 3D CO study of 260 ETGs (Young 
et al. 2011), though the ATLAS 3D sample has a larger 
volume out to distances 46 Mpc (compared to 25 Mpc for 
the HRS). Comparing those HRS galaxies which are in 
both the ATLAS 3D and HRS (i.e., observed in both CO 
and dust continuum), the detection fraction in CO is a 

22 Using the two-sided Fisher's exact test, the equivalent of the 
X-squared test but for small samples. 



factor of two lower than at 250 fim. Similarly for those 
galaxies with Herschel and Hi observations, the detec- 
tion fraction in Hi is ~1.5 times lower (Table|3]). The 
Herschel observations of the continuum emission from 
the dust are therefore currently the most sensitive way 
of detecting the ISM in the HRS sample of ETGs. 

Given the small number of ellipticals that have been 
detected by Herschel, it is interesting to ask the question 
of what makes these special. Twenty of the ETGs in 
our sample are classified as having low-ionization narrow 
emission line regions (LINERs) . The origin of LINER ac- 
tivity in ETGs can be ionizing photons from the old stel- 
lar population (e.g., di Serego Alighieri et al. 1990; Sarzi 
et al. 2010), heating by the hot ISM (Sparks et al. 1989) 
or AGN heating (e.g., Gonzalez-Martin et al. 2009). At 
least 12 of the HRS sample have bright radio cores (Fig- 
ure [1} suggesting AGN photoionization is important for 
at least 20% of the sample. For the most massive el- 
lipticals, the hot ISM could be a significant photoioniz- 
ing source; nine of the detected LINER sources in our 
sample have X-ray luminosities well above that expected 
from the discrete stellar population i.e., Lx/Lb ~ 10 29 5 
(O'Sullivan et al. 2001; Section 14. 2\ suggesting that a 
significant fraction of the LINER emission may originate 
from the hot ISM rather than an AGN source. Of the 
seven detected ellipticals, two are classified as LINERs, 
two as Seyferts and a further two are unambiguously 
identified as containing a bright AGN (Ho et al. 1997; 
Schmitt 2001; Veron-Cetty & Veron 2010). Only one of 
the seven, NGC 4406 (HRS150) is a 'pure' elliptical with 
no sign of an active nucleus, but this galaxy has acquired 
dust and atomic gas through a recent interaction (Gomez 
et al. 2010). Nine other 'pure' ellipticals in our sample 
are not detected at 250 /xm. 

Another way to look at this question is to compare 
the properties of the detected and undetected sources in 
other wavebands. The distribution of the £>-band lumi- 
nosities for the sample is shown in Figure|3]along with the 
distributions of the detected sources. The figure shows 
that while there is no tendency for the detected SO galax- 
ies to be optically luminous, the detected elliptical galax- 
ies are optically the most luminous ones. 

Figure[3]shows a similar comparison for the X-ray lumi- 
nosities. In this case, one might suspect that the galax- 
ies with the highest X-ray luminosities from the hot halo 
are the least likely to be detected by Herschel because of 
sputtering. The figure shows that the elliptical galaxies 
detected at 250 /im tend to have higher X-ray luminosi- 
ties so that elliptical galaxies more luminous in X-rays 
and optical (i.e., the most massive stellar systems) tend 
to be detected by Herschel. The detected SO galaxies 
span a wide range of Lb and Lk, indicating that the 
stellar mass of these systems is not the critical factor for 
detection by Herschel. 

In contrast to the X-rays, we find that galaxies con- 
taining large masses of cool gas are more likely to be 
detected by Herschel (Figure [3]). CO observations of our 
sources taken from the ATLAS 3D survey (Young et al. 
2011) suggest that 28% of our ETGs have molecular gas 
(Table 131, similar to the detection rate of 22% for the 
entire ATLAS 3D sample of 260 ETGs. Atomic hydrogen 
observations are available in the literature for 79% of our 
sample, although only 35% of these galaxies are detected. 
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Table 3 

The Number of Early-type and Elliptical Galaxies Detected with Herschel 



Type SPIRE Sample X-Ray Sample CO Sample 



Observed Detected 250 Observed Detected X-ray Detected 250 Observed Detected CO Detected 250 

N % N % N % N % N % N % N % N % 



Total 


62 


100% 


31 


50% 


38 


61% 


29 


76% 


15 


52% 


56 


90% 


16 


28% 


16 


100% 


E 


21 


34% 


5 


24% 


20 


53% 


16 


80% 


4 


25% 


20 


34% 





0% 





0% 


SO 


38 


61% 


23 


63% 


15 


39% 


11 


77% 


9 


82% 


35 


63% 


16 


55% 


16 


100% 


pec 


3 


5% 


3 


100% 


3 


8% 


2 


67% 


2 


100% 


1 


2% 





0% 





0% 



Hi sample Environment 



Observed Detected Hi Detected 250 Obs. Virgo Det. Virgo Obs. Not Virgo Det. Not Virgo 
N % N % N % N % N % N % N % 



Total 


49 


79% 


17 


35% 


16 


94% 


47 


76% 


21 


45% 


15 


24% 


10 


66% 


E 


15 


31% 


1 


27% 


4 


100% 


17 


74% 


5 


29% 


6 


26% 


2 


33% 


SO 


31 


69% 


13 


38% 


12 


92% 


30 


77% 


16 


53% 


9 


23% 


8 


88% 


pec 





0% 





0% 





0% 



















(b/a < 0.5) (b/a > 0.5) 

Observed Detected Det. Virgo Det. Not Virgo Observed Detected Det. Virgo Det. Not Virgo 

N % N % N % N % N % N % N % N % 



SO 18 46% 12 67% 13/17 76% 1/1 100% 21 54% 12 57% 2/10 20% 9/11 82% 

Barred 11 28% 6 55% 5/10 50% 1/1 100% 13 33% 6 46% 4/10 40% 2/3 67% 

Notes. — Top Panel: the detection numbers for the 250 /^m, X-ray, and CO samples in this work. Columns arc as follows. Column 1: morphological 
type for the sample. Columns 2— 5:SPIRE-c\ctcctcd sources including — Column 2: number of sources observed; Column 3: percentage of the observed 
sample which falls into the different types; Column 4: number of sources detected with SPIRE; Column 5: percentage detected with SPIRE compared to 
the observed sources in each morphological class. Columns 6— ll:X-ray sources including — Columns 6 and 7: number of sources in our sample observed 
in X-ray from the literature and the percentage which falls into each morphological class (C Sullivan & Ponman 2004; Pellegrini 2010); Columns 8 and 
9: number and percentage of the observed X-ray sources which are detected in X-ray; Columns 10 and 11: number of sources detected in X-rays and 
detected with SPIRE. Columns 12-17: CO sources including — Columns 12 and 13: number of sources in our sample observed in CO from the literature 
and the percentage which falls into each morphological class (Young ct al. 2011); Columns 14 and 15: number and percentage of the observed CO sources 
which are detected in CO; Columns 16 and 17: number of sources detected in CO and detected with SPIRE. Middle Panel: the detection numbers for 
the Hi and 250 nm sample, and for the entire 250 fim sample split by environment. Columns arc as follows. Columns 2-7:Hl sources including — Columns 
2 and 3: number of sources in our sample observed in Hi from the literature and the percentage which falls into each morphological class; Columns 4 
and 5: number and percentage of the observed Hi sources which arc detected in Hi; Columns 6 and 7: number of sources detected in Hi and detected 
with SPIRE. Columns 8-15: Environment: the number and percentage in each subset defined as within the Virgo Cluster (Columns 8—11) and outside 
(Columns 12-15). Given the small numbers in the Virgo/non- Virgo samples, we do not separate the peculiar galaxies in this analysis. Bottom Panel: 
the detection numbers for SO+SOa and barred SO galaxies split into regimes depending on their axial ratios, where b/a < 0.5 and b/a > 0.5, indicating 
edge-on or face-on galaxies, respectively. Columns are as follows. Columns 2-9:Edge-on galaxies including — Columns 2 and 3: number and percentage 
of edge-on sources (b/a < 0.5); Columns 4 and 5: number and percentage of these sources which are detected with SPIRE. Columns 6 and 7: number 
of sources detected with SPIRE inside the Virgo Cluster and Columns 8 and 9 for those outside Virgo. Columns 10— 17:-Face-on galaxies — Columns 10 
and 11 - number and percentage of face-on sources (b/a > 0.5): Columns 12 and 13: number and percentage of those sources which arc detected with 
SPIRE. Columns 14 and 15: number of sources detected with SPIRE inside the Virgo Cluster and Column 16 and 17 for those outside Virgo. 



We can also consider whether certain morphological 
types are more likely to be detected. Eighty percent of 
the undetected galaxies are morphologically classified as 
barred systems (NED, Table [I]), which at first glance 
suggests that barred galaxies have less dust content than 
galaxies without a bar. However, we actually detect 50% 
of the barred sources from the HRS sample (Table [3]), 
suggesting that the presence of a bar does not determine 
whether it is detected at 250 /im. Finally, we investigated 
whether the inclination angle might affect whether or not 
the galaxy was detected by Herschel using the measured 
axial ratios of the galaxies (minor axis/major axis — b/a). 
The detection rate for edge-on galaxies (b/a < 0.5) is 
67% and for face-on galaxies (b/a > 0.5) is 57% though 
this is not statistically significant. Whether a galaxy has 
rings, bars or high inclination angles therefore does not 
affect whether it is detected by Herschel. 

One important effect on whether we detect ETGs with 
Herschel is environment. It has long been known that 
galaxies in dense environments are depleted in atomic 
gas, and recent Herschel observations have shown that 
this is also true of dust (Cortese et al. 2010b; C12). 
Gravitational interactions might also either remove or 



add gas and dust to a galaxy. We can do a rough com- 
parison of the detection rates in different environments 
using the cluster or cloud membership listed in Table [TJ 
Following C12, we split our sample into two sets: one 
set includes galaxies inside the Virgo Cluster and its 
outskirts (N = 47 galaxies) and the other set includes 
galaxies outside Virgo (N = 15). We detect 53% (16/30) 
of the SOs and 29% (5/17) of the ellipticals within the 
Virgo Cluster. Outside Virgo, we detect 89% (8/9) and 
33% (2/6) respectively (Tabled]). The detection rates 
are therefore lower for galaxies within the Virgo Cluster, 
but this is not at a significant level ™ 

3.2. Residual Star Formation? 

We now discuss the NUV-r color for our sample; for 
star-forming galaxies, this is a useful measure of the spe- 
cific star formation rate (SFR/A/*) since the NUV traces 
recent star formation and the r band is a proxy for stel- 
lar mass (e.g., Wyder et al. 2007). For Milky Way like 
extinction curves, the NUV is less affected by internal 
dust obscuration compared to the FUV, but the NUV 
can also be affected by the presence of evolved hot stars 
(the UV-upturn, e.g.. Greggio & Renzini 1990; Yi et 
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Figure 3. Top left: the range of B-band luminosities; top right: X-ray luminosities for the SO+SOa (top) and elliptical (bottom) galaxies. 
The outer dashed histogram shows the range of the whole sample (SO+SOa or E combined), the thick black unfilled histograms are the 
subsamples (SO+SOa and E separated out), and the detected sources are shaded gray. Bottom left: the molecular gas mass estimated from 
CO observations; bottom right: the atomic gas mass from Hi. The red histograms show the upper limits on the gas mass obtained from 
the CO (Hi), i.e., these are sources in our sample not detected either at 250 ^m or CO (Hi). The red shaded histogram shows the sources 
detected at 250 fim but not detected in CO (Hi). 



al. 2011). Blue, star-forming galaxies occupy the re- 
gion NUV-r < 3.5 (Kaviraj et al. 2007), the so-called 
'blue cloud', with the red, passive ETGs found at NUV- 
r > 4.5 in the "red sequence". Although the color cut 
of NUV-r > 4.5 is often used to select galaxies on the 
red sequence, the most passive galaxies are significantly 
redder than this. Saintonge et al. (2011) have found that 
galaxies with NUV-r > 5 have very little molecular gas 
and Kaviraj et al. (2007), Schawinski et al. (2007), and 
Yi et al. (2005) have showed that quiescent non-star- 



forming ETGs occupy the region above the boundary 
NUV-r = 5.4. 

Figure [4] shows the NUV-r color plotted against stel- 
lar mass for all the HRS galaxies taken from C12 and 
L. Cortese et al. (in preparation). Sources detected at 
250 (im are highlighted with the outlined black circles. 
All but two of our early-type sample have NUV-r > 4.5 
and most cluster around the line NUV-r = 5.5. If 
NUV-r < 5.4, the blue colors may be due to recent 
star formation but they may also be caused by the UV 
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upturn. NGC4552 is a famous example of this, with 
NUV-r ~ 5.2, even though there is no residual star 
formation (Figure HJ). Excluding NGC4552, 12 out of 
18 SOs which lie below this boundary are detected with 
Herschel, and 2 out of 7 ellipticals. The two ellipticals 
are the AGN NGC4486 (HRS138) and the (unusual) X- 
ray source NGC4636 (HRS241); the NUV emission from 
these galaxies may be related to their central activity. 

As we detect dust emission in almost half of the galax- 
ies, the UV fluxes should be corrected for the dust at- 
tenuation since this will cause the galaxies to be redder 
in NUV-r than they truly are. Given the complexity of 
disentangling the old and young stellar populations, we 
leave this for future work and, given the sensitivity of the 
NUV emission to internal dust extinction, we make no 
attempt to estimate the SFRs from the NUV-r colors. 
Following Schawinski et al. (2007), we simply classify 
galaxies as quiescent where NUV-r > 5.4 and possibly 
undergoing recent star formation (RSF) where NUV- 
r < 5.4. In this case (including sources not detected 
with GALEX but excluding NGC4486 and NGC4552), 
~ 40% of our early-type sample is in the RSF group, 
including 30% of our ellipticals. This agrees with the 
fraction of SDSS ETGs found with RSF (Schawinski et 
al. 2007). Of course, our estimate may be biased in a 
number of ways, since we have not corrected for dust ex- 
tinction and the colors of some of the galaxies in the RSF 
sample may be due to LINER-type activity or to the UV 
upturn. 

The most interesting result, however, that can be seen 
in Figure 2] is the very large number of ETGs that lie on 
the red sequence but are still detected by Herschel. It 
is often assumed that the red sequence contains galaxies 
with virtually no current star formation and little inter- 
stellar gas, but the detection with Herschel of ~ 50% of 
these galaxies clearly shows that they contain significant 
reservoirs of interstellar material agreeing with recent ob- 
servations of CO in ETGs (Young et al. 2011). 

3.3. Dust Masses and Temperatures 

The global IR-submm SEDs of the detected 
S0+S0a and elliptical galaxies with available data 
from 60 to 500 ^m are presented in Figures [5] and [6l 
respectively. NGC4486 (M87) and NGC4374 are both 
bright radio sources and their SEDs are well fitted by 
a power law at the longer FIR wavelengths. For the 
former, we find no evidence for dust emission above the 
power-law synchrotron component (in agreement with 
Baes et al. 2010). In NGC4374, we do see a strong 
excess from dust emission in the FIR. We fit the SEDs 
from 100 to 500 with a modified blackbody model, 
where 

_ Kv M d B(y,T d ) 

M d is the dust mass, T d is the dust temperature, B(u, T d ) 
is the Planck function and D is the distance to the galaxy. 
K v is the dust absorption coefficient described by a power 
law with dust emissivity index /3, such that jc„ oc v" . For 
/3 = 2 (typical of Galactic interstellar dust grains), we 
use K35o M m = 0.19m 2 kg -1 (Draine 2003). Although k is 
notoriously uncertain, virtually all of our analysis relies 
on k being constant between galaxies rather than on its 
absolute value. 



We assumed D = 17Mpc for the Virgo Cluster 
(Gavazzi et al. 1999) and 23Mpc for the Virgo B cloud; 
for all other galaxies, we estimated D from the re- 
cessional velocities using a Hubble constant of H a = 
70 kms _1 Mpc _1 . 

We found the best-fit solution by minimizing the chi- 
squared (x 2 ) difference between the model and the mea- 
sured fluxes, involving the model with the filter trans- 
mission functions. To account for the uncertainties in 
the fluxes that are correlated between the SPIRE bands 
(5%), we used the full covariance matrix in the x 2 cal- 
culation. The covariance matrix has diagonal elements 
with the total variance of each band and non-diagonal el- 
ements consisting of the covaraince between bands from 
the correlated calibration uncertainties. We treated flux 
measurements at wavelengths < 70 //m as upper limits 
for the cold dust emission as previous works (Smith et 
al. 2010; Bendo et al. 2010a) found that there is a signif- 
icant contribution at these wavelengths from a warmer 
component of dust. We estimated uncertainties in the 
dust temperatures and masses using a bootstrap tech- 
nique, in which we generated artificial sets of fluxes from 
the measured fluxes and errors, and then applied our fit- 
ting technique to each set of fluxes. The dust tempera- 
tures, and therefore dust masses, depend on the choice of 
/3; = 1.5 would yield slightly higher dust temperatures 
for the whole sample (Bendo et al. 2003) and therefore 
lower dust masses. 

The SED for NGC5273 (Figure E} is the only source 
which appears to have an additional component present. 
Using the (3 = 2 one-temperature component described 
above provides an adequate fit to the photometry data 
if an additional synchrotron component is present. The 
SED is also adequately fit with a one temperature dust 
component model with (3 as a free parameter, (3 = 0.9, 
or with a two-temperature dust component model. To 
fit the data with the latter model would require an addi- 
tional component of ~ 9 K dust which would increase the 
dust mass by an order of magnitude. Both the (3 = 2 + 
synchrotron model and (3 = 1 model give a dust mass of 
10 5 45 Mq which is somewhat lower than the dust mass 
estimated using the (3 = 2 one-temperature dust com- 
ponent without synchrotron (10 5 65 Mq) but within the 
errors (Table S|) . As no millimeter data are available for 
this source, we cannot rule out synchrotron contamina- 
tion, so in this analysis we have chosen to stick with the 
(3 = 2 one-component model to be consistent with the 
other sources. 

We have estimated upper limits to the dust masses for 
the galaxies undetected by Herschel using the flux upper 
limit at 250 fim (including the noise on the image, an 
estimate of confusion noise and background variance in 
the map as described in L. Ciesla et al., in preparation) 
and the average temperature for the detected galaxies. 

The dust masses and temperatures are presented in 
Table [H and we show the range of dust temperatures, 
masses, and stellar masses for the SO+SOa and elliptical 
detected sample in Figure [7] These range from M d = 
10 5.o-7.i Mq and ie_32 K with mean values of 10 614 M 
and 24 K, respectively (Table [5]). The stellar masses of 
these galaxies (C12) range from io 9 - 9-11 - 2 M & with an 
average of (logM*(M©)) = 10.89; the average dust-to- 
stellar-mass ratio is \og(M d /M r ) = -4.33 ± 0.14. The 
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Figure 4. Left: the NUV-r color vs. the stellar mass of the HRS sample (L. Cortese et al. in preparation) including our morphologically 
classified spirals (gray, C12), SO+SOa (blue), and elliptical galaxies (red). The SO+SOa and elliptical galaxies detected with SPIRE are 
highlighted with the outer black circles. Note that the galaxy NGC4486 is shown as detected in this plot but we do not see evidence for 
dust (above the synchrotron emission), NGC4486 is highlighted with "+". The well-known galaxy with no evidence for star formation 
(NGC4552) but with somewhat bluer NUV-r color than the average elliptical in our sample is highlighted with a "x". UV photometry 
exists for 87% of the ellipticals and 95% of the SO sample. Right: the distribution of NUV— r colors for the observed samples split into 
spirals (gray), SO+SOa (blue), and ellipticals (red; unfilled histograms) including detections for the early -types (shaded histograms). 



Md /M* ratios are plotted in Figure |5] and clearly shows 
that, at the same stellar mass, there is a sharp fall in the 
dust-to-stellar-mass as we move from spirals to SOs (by 
roughly a factor of 10), and that this fall continues as we 
move from SOs to ellipticals (again, by a factor of 10). 
The two ellipticals with anomalously high values given 
their stellar mass are NGC3226 (HRS3) and NGC4406 
(HRS 150). Both galaxies show signs of a tidal interaction 
(see Appendix , suggesting that the dust may have 
been acquired as the result of an interaction with a dust- 
rich galaxy. 

4. DUST, STARS AND GAS 

In this section we compare the dust temperatures, 
masses, and dust-to-stellar ratios with other Herschel re- 
sults of late-type and early-type galaxies. First, the tem- 
peratures for the HRS ETGs are generally higher than 
temperatures for the galaxies detected in blind submm 
surveys or in other samples of nearby galaxies (see Planck 
collaboration 2011; Dunne et al. 2011) when corrected for 
different choices of (3. We have tested whether the tem- 
peratures of ETGs are systematically higher than those 
of late types by comparing our dust temperatures with 
those of 71 Virgo galaxies in Davies et al. (2012) who 
used an identical method. We find that the tempera- 
tures of the early types are systematically higher (with 
U = 1638, m = 70, n 2 = 30 and P > 99.999Q, in 
agreement with the result that the ETGs in the HRS have 
warmer IR colors than late types (Boselli et al. 2010a; 
see also Bendo et al. 2003; Engelbracht et al. 2010). 

In Figure we compare the ratio of dust-to-stellar- 
mass for the HRS ETGs with the Herschel KINGFISH 

23 Using the Mann— Whitney U statistic, appropriate when com- 
paring two groups where the underlying distribution of the data is 
not necessarily normal. The statistic assumes that the observations 
are independent and are continuous (i.e., able to be ranked) and is 
more robust than the Student's t-test. 



(Kennicutt et al. 2011) results based on 10 nearby early 
types (Skibba et al. 2011). They find warmer dust tem- 
peratures, with a mean of (T^) = 30 K (mostly due to 
the use of a lower emissivity index ft = 1.5 compared to 
this work), and their sample has a lower average stellar 
mass compared to the HRS. Though their dust masses 
are similar, the dust-to-stellar-mass ratios are higher in 
KINGFISH than for the HRS. This may be a result of 
the different environments (a large fraction of the HRS 
galaxies are in a rich cluster) and/or the selection of the 
sample (KINGFISH ETGs may include unusual galaxies 
rather than the flux-selected ETGs in the HRS). 

Rowlands et al. (2012) detected only ~ 5% of ETGs in 
a blind submm survey with Herschel at 250 /im, finding 
a median dust mass for their detected galaxies (in their 
lowest redshift bin) of 10 78 M© with average dust tem- 
perature consistent with the HRS ETG sample. Their 
dust mass is significantly higher than the dust masses 
of the HRS ETGs and their mean value of their dust- 
to-stellar mass ratio ( log (Md/M*) = —2.95, Figure [Sj> is 
also larger. Many of their sources have bluer UV-optical 
colors than the HRS (with a significant fraction lying 
below NUV-r < 4: their Figure 15), suggesting that the 
larger amount of dust in the H- ATLAS galaxies is as- 
sociated with increased star formation (Rowlands et al. 
2012). The simplest explanation of the differences is that 
the H-ATLAS study is picking up the rare, very dusty 
ETGs, which the HRS misses, partly because it is not 
large enough, and possibly also because it is dominated 
by the ETGs in the Virgo Cluster. A possible local ex- 
ample of these rare dusty ETGs is the local elliptical Cen 
A (T. Parkin et al. 2011, submitted). Since the depth of 
the SDSS images used for the morphological classification 
in the H-ATLAS study made it impossible to distinguish 
between ellipticals and SOs, it is also possible that some 
of the ETGs detected in H-ATLAS are SOs or early-type 
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Figure 5. Spectral energy distributions of the 24 detected SO+SOa galaxies in our sample. The red solid line is the best-fit modified 
blackbody, blue points are fluxes used as upper limits, and black crosses are the photometry used for the cold dust component, including 
Herschel SPIRE fluxes and PACS, IRAS, and Spitzer where available. Errors shown also include calibration errors (L. Ciesla et al., in 
preparation). The best-fit model parameters are provided in Table [4] 



spirals rather than ellipticals. 

Rowlands et al. (2012) were also able to estimate the 
dust masses of the ETG population as a whole using a 
stacking analysis. They estimated a mean dust mass of 
10 6 - 3 M Q , with a mean dust temperature of 25 K and a 
mean dust-to-stellar-mass ratio of log(Md/M*) = —4.87. 
This is in reasonable agreement with the average param- 
eters of the early-type sample from this work and in par- 
ticular, with the SO+SOa sample and our most dusty el- 
lipticals (Figure [5]). 

4.1. Dust and stellar mass along the Hubble Sequence 

The average dust mass for the SO+SOa and elliptical 
galaxies detected by Herschel is 10 63 M Q and IO 5 5 M© 



(Figure Table [5]) with average dust temperatures of 
23.5 and 25.7 K, respectively. The dust masses for the 
two morphological groups are significantly different ( 
Mann-Whitney £7=117 with n x = 6, n 2 = 24 at P > 
99.7% for two-sided test)P^ This does not take into ac- 
count the large number of upper limits so we have used 
the Astronomical Survival Analysis programs (Feigelson 
& Nelson 1985), implemented through the IRAF STS- 
DAS statistics package, to compare the two samples. 
Using three different tests we find that the probabilities 
that the elliptical and SO+SOa dust masses are drawn 
from the same population are 0.0005, 0.0006, and 0.004, 
and for M^/M*, the probabilities are zero (Table [5]); 
these parameters for SO+SOa and ellipticals in our sample 
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Figure 6. Spectral energy distributions of the seven detected ellipticals in our sample. Blue points are those fluxes used as upper limits 
and black crosses are the photometry used for the cold dust component, including Herschel SPIRE and PACS, IRAS, ISO, and Spitzer 
where available. Radio fluxes are also shown by black points for NGC4374 and NGC4486 (M87). The red solid line is the best fit to the 
data and includes a modified (single temperature component) blackbody with (5 = 2 and a synchrotron component. Errors shown also 
include calibration errors (L. Ciesla et al. in preparation). These galaxies are discussed in more detail in Appendix |A1 and the best-fit 
model parameters are listed in Table [4] 

almost certainly arise from different distributions. We 
have used the Kaplan-Meier estimate^ to estimate the 
mearpl for the two samples, making the necessary ap- 
proximation that the lowest upper limit is actually a de- 
tection. With this approximation, we find that the mean 
dust mass for the SOs is 10 5 ' 87 M Q and the mean dust 
mass for the ellipticals is 10 5 21 M Q (Table©. 

C12 has recently shown that the dust-to-stellar-mass 
ratio decreases when moving from late- to early-type 
galaxies, but used a conservative approach by treat- 
ing non-detections as upper limits. Here we incorpo- 
rate the information in the upper limits on the non- 
detected sources and quantify the change as we move 



24 It is likely that the censoring in this data set is not random: 
the censored data are a multiple of the noise in the data and we 
have observed to a flux limit in K with shorter integration times 
for the SO+SOa sample compared to the ellipticals. The censored 
data may have some randomization introduced by the SED-fitting 
technique since the SED can have different distributions but still 
give the same M$. We note that the censored data could be biased. 

25 The median of the distribution function returned by Kaplan- 
Meier estimator is well defined if the censored data are restricted 
to less than half of the sample. Given the small sample size of the 
ellipticals which is dominated by upper limits (70% of the data) 
we use the mean in this work. 



along the Hubble sequence. We estimate that the mean 
value^ of the dust-to-stellar-mass for the ellipticals is 
log(M d /M*) = -5.83 ±0.1 and for the SOs is -4.42 ±0.1 
(Table [5j. 

Following C12 (his Figure 5), we plot the dust-to- 
stellar-mass ratio against morphological type in Fig- 
ure [S] We split the late- type sample into the subgroups: 
Sa±Sab, Sb+Sbc, Sc+Scd and Sd±Sdm. The figure also 
shows the strong drop in mean dust-to-stellar mass ratio 
from early-type spirals to SOs (and the further decline 
as we move to ellipticals) . Figure |9] suggests that the 
variation in this ratio is larger for early types than late 
types. This is particularly true for early-type spirals and 
S0±S0a galaxies with Md/M* ranging by a factor of ap- 
prox. 50. Although ellipticals probably have a similarly 
wide range in dust-to-stellar-mass, the high number of 
non-detections make it difficult to be sure of this. 

4.2. The origin of dust in ETGs 

One of the big questions about dust in ETGs is its ori- 
gin. There are three main possible sources: (1) the dust 
is formed in the old stellar population in the atmospheres 
of AGB stars, (2) the dust is acquired externally as the 
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Notes. — Column 1: HRS ID; Column 2: morphological classification; Columns 3 and 4: the 
FIR luminosity (of the cold dust component) and dust temperature; Columns 5: dust masses 
estimated using j3 — 2 and K350 — 0.19m 2 kg -1 , 3 a upper limits (L. Ciesla et al., in prepara- 
tion) arc quoted where the source is not detected using the mean dust temperature; Columns 
6 and 7: molecular masses estimated from CO (Young et al. 2011) and HI masses (references 
below); Columns 8: slow (S) or fast (F) rotator as defined by ATLAS 3D (Emscllcm et al. 2011). 
References. 

1 Gavazzi ct al. (2003) 

2 Haynes ct al. (2011). 

3 Springob ct al. (2005). 

4 Noordcrmccr ct al. (2005). 

a No SDSS data arc available for HRS209 or HRS186. 
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Figure 7. Top: histogram of the dust temperatures based on sin- 
gle modified blackbody fits; bottom: dust (left) and stellar masses 
(right shaded; L. Cortese et al., in preparation) for the detected 
ellipticals (red) and SO+SOa galaxies (blue) — see Table|4] The pop- 
ulations of the dust and stellar masses of the detected ellipticals 
and SOs arc different with significance level P > 97.9% (Mjj) and 
P > 99.9% (M») using a Mann- Whitney U-test. 

result of a galaxy merger or a tidal interaction, and (3) 
the dust has the same rather uncertain sources as in late- 
type galaxies (e.g., Morgan & Edmunds 2003). In (3), 
the dust might be continuously created by grain growth 
in the ISM or current star formation via dust production 
in supernovae (e.g., Barlow et al. 2010; Matsuura et al. 
2011; Gomez et al. 2012), or it might be left over from a 
more vigorous star-forming epoch. The strongest predic- 
tion is made by the first hypothesis, since this makes the 
clear prediction that there should be a strong correlation 
between dust mass and stellar mass. 

To test point (1) we compare the FIR emission from 
the dust with the starlight and the emission from the hot 
ISM. If the source of the hot gas is mass loss from stars, 
one would expect the X-ray luminosity (Lx) to correlate 
with the optical luminosity (Lb) which is true for our 
sample (Appendix IB1 Equation (|B1[) ) ■ Similarly, if the 
source of the dust is stellar mass loss, we would expect 
to see a correlation between dust mass and optical lumi- 
nosity; alternatively, a lack of correlation may indicate 
that mergers are important for delivering dust. 

The FIR luminosity, Lfir (Table 2]) was determined 
by integrating over the SED (Figures [5] and [6]) . (This is 
an underestimate of the total IR luminosity as the 70 /im 



fluxes are treated as upper limits in our SED fitting; 
we estimate that including a warm temperature com- 
ponent to fit the MIR-70 /jm emission would contribute 
(on average) 3% to the FIR luminosity, though for some 
sources this can be up to 17% — see also Muhoz-Mateos 
et al. 2009.) The mean ratio of FIR to optical luminos- 
ity log(L Fm /L B ) for SOs is -1.1 ±0.1 and -2.2 ±0.1 for 
ellipticals, therefore ellipticals have less FIR luminosity 
per unit blue luminosity than SOs (also seen in an IRAS 
sample of isolated galaxies - Lisenfeld et al. 2007). In 
Figure I10[ we compare the -Lfir with the optical lumi- 
nosity which does indeed increase as Lb increases, but 
the large scatter suggests that dust may have been ac- 
quired from tidal interactions or mergers as well as stellar 
mass loss. There is no evidence of a correlation between 
Lfir and Lb (Figure [TUl Appendix iBj) . We also find no 
correlation with Lfir and the X-ray luminosity Lx ■ 

Since Lfir depends strongly on the temperature of 
the dust, it is instructive to compare the dust mass, 
Ma, with t he o ptical and X-ray luminosities (Figure [TUl 
Equations (|B4[) and (|B5[)). We find that the dust mass 
increases weakly with Lb but there is no statistical cor- 
relation between the two. The same is true for Md 
and Lx- We compare the predicted relationships for 
a model in which dust is replenished by stellar mass 
loss from evolved stars and destroyed by dust sput- 
tering (Goudfrooij & de Jong 1995; Figure [TOj). The 
mass of gas lost in these winds is estimated by M gas = 
1.5 x lO- n (L B /L ) A-foyr" 1 (Faber & Gallagher 1976), 
which ranges from 0.005 to 0.75 Mgyr" 1 for the HRS 
ellipticalfO- The estimates of the final dust masses de- 
pend critically on the assumption of the lifetime of the 
dust: dust destruction in a hot plasma probably occurs 
on timescales of Td ~ 10 6 -10 7 yr (Draine & Salpeter 
1979). Goudfrooij & de Jong note that Td ~ 10' yr is 
the maximum plausible destruction timescale since this 
corresponds to electron densities in the hot X-ray plasma 
«h < 10~ 3 cm~ 3 . In the figure we plot the predicted re- 
lationships for values of log 10 (r ( j) from 6.5 to 7.5. The 
ETG sample detected by Herschel have dust masses well 
above the theoretical curves, suggesting an alternative 
source of dust is needed for most of the sources. Longer 
grain lifetimes for some of the galaxies may also occur 
if a significant component of cool gas exists (see Section 
14. 3|k in this case dust may reside in conditions similar 
to the Milky Way, with destruction timescales of the or- 
der of log 10 (Td) = 8.5. Three elliptical galaxies have dust 
masses consistent with the boundary marked by the max- 
imum destruction timescale log 10 (r ( i) = 7.5. Note that 
the figure shows that dust originating from AGB-stellar 
mass loss (Clemens et al. 2010) could be consistent with 
the upper limits estimated for 50% of the ETGs. 

We conclude then that there is no strong statistical 
evidence of a correlation in FIR luminosity or dust mass 
with blue- or X-ray luminosity. The dust masses of the 
ETGs detected by Herschel are also larger than the es- 
timates from a model in which dust is produced from 
stellar mass loss from evolved stars, even with the most 
generous assumption about how long the dust grains will 
survive in the ISM. Furthermore, if the dust in ETGs is 

26 We assume that the gas-to-dust ratio in the stellar winds of 
evolved stars is ~ 150. 
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Figure 8. Dust-to-stellar-mass ratio vs. stellar mass. Blue points are the detected SO+SOa galaxies, with red points for the ellipticals. 
Upper limits on the dust mass are shown with triangles. The mean M^/M* for the HRS spirals in each stellar mass bin is shown with 
gray circles (see C12, his Figure 5) with error bars indicating the error on the mean. We also compare the mean M^/M* results of late 
types in HeViCS (blue star; Davies et al. 2011) and H-ATLAS (yellow star; Dunne et al. 2011), early types with H-ATLAS (Rowlands et 
al. 2012 for detected (purple diamond) and undetected (cyan diamond) early- type populations), and early types with KINGFISH (purple 
star; Skibba et al. 2011). The stellar masses are from L. Cortese et al. in preparation. 



Table 5 

Mean Parameters for the Sample 



Type 


N 


T d 
(K) 


logM d 

(M @ ) 


logM* 

(Me) 


logM d /M, 


logLx 
(ergs" 1 ) 


logAf H2 
(M e ) 


logAfni 

(Me) 


Detected sample: 
Total 

E (+E/S0+pec) 
SO (+pec) 


31 
7 
24 


23.9 ± 0.8 
25.7 ± 2.5 
23.5 ± 0.8 


6.12 ± 0.12 
5.54 ± 0.28 
6.26 ± 0.13 


10.49 ± 0.07 
10.89 ± 0.19 
10.40 ± 0.06 


-4.33 ± 0.13 
-5.26 ± 0.24 
-4.13 ± 0.12 


40.02 ± 1.17 (15) 
40.96 ± 3.11 (5) 
39.34 ± 1.90 (10) 


< 7.21 (0) 
7.84 ± 0.25 (16) 


8.40 ± 0.17 (15) 
8.26 ± 0.14 (5) 
8.47 ± 0.19 (10) 


Entire sample: 
Total 

E (+E/S0+pec) 
SO (+pec) 


62 
22 
39 




5.59 ± 0.09 
5.21 ± 0.09 
5.87 ± 0.13 


10.09 ± 0.12 
10.66 ± 0.13 
10.38 ± 0.04 


-5.12 ± 0.14 
-5.83 ± 0.11 
-4.42 ± 0.10 


39.89 ± 0.18 
40.13 ± 0.24 
39.25 ± 0.24 


7.13 ± 0.67 (57) 
7.02 ± 1.29 (21) 
7.25 ± 0.75 (36) 


7.57 ± 0.11 (48) 
7.72 ± 0.17 (15) 

7.58 ± 0.13 (33) 



Notes. — Top: the mean parameters for the detected sample of early types. Bottom: the mean parameters estimated using survival analysis, 
including the upper limits. The error quoted is the standard error of the mean. The columns are as follows. 

Column 1: morphological type. Column 2: the number of sources. Column 3: the average dust temperature. Column 4: the average dust mass 
(dust masses are estimated using a single temperature modified blackbody with — 2 and K350 — 0.19m 2 kg~ 1 ). Column 5: the average stellar 
masses using the optical colors (L. Cortese et al., in preparation). Column 6: the average dust-to-stellar-mass ratio. Column 7: average X-ray 
luminosity — the numbers in brackets give the number of sources with X-ray and 250 /xm detections. Column 8: the average molecular gas masses 
estimated from the CO data. Column 9: the average atomic gas masses estimated from Hi data. 



in equilibrium, with dust formation in old AGB stars 
balanced by dust destruction via sputtering, we would 
also expect a small range in M^/M*. Figures [8] and [9] 
show that the variation in this ratio is much larger for 
early types than late types; this is particularly true for SO 
galaxies with M^/M* ranging by a factor of ~50. Finally, 
another argument against hypothesis (1) and the super- 
nova origin of dust in hypothesis (3) is that if stellar mass 
loss was responsible for the dust, the dust would be dis- 
tributed in a similar way to the stars. However, for five 
of the six ellipticals in which Herschel has detected dust 
(excluding M87), there is evidence from the appearance 
of the dust in absorption against the optical or near-IR 
continuum emission that the dust is distributed in a dif- 



ferent way from the stars (see Appendix [SJ. 

The evidence therefore points toward a merger origin 
of dust in this sample, this is supported by evidence from 
other recent studies that the cool ISM in ETGs may have 
been acquired as the result of gravitational interactions. 
For example, Davis et al. (2011) have used the misalign- 
ment between the kinematics of the stars and the gas in 
ETGs to conclude that at least one-third of ETGs have 
acquired much of their ISM by this means. 

4.3. A cool interstellar medium in early-type galaxies? 

In the last decade, it has become clear that a significant 
fraction of ETGs do contain a cool ISM similar to that 
of spirals. Morganti et al. (2006) detected Hi in 70% of 
their sample of nearby ETGs (see also Oosterloo et al. 
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E SO+SOa Sa+Sab Sb+Sbc Sc+Scd Sd+Sdm 

Hubble Type 

Figure 9. Dust-to-stcllar-mass ratio vs. morphological type for the detected and undetected sources. Gray circles show the M<j/Af* 
estimated for detected sources, with upper limits highlighted with gray triangles. The mean and standard deviation of M^/M* for the 
detected samples are highlighted with black boxes and error bars. Offset from the detected sources, we also show the entire HRS sample in 
each Hubble type including both detected sources and upper limits. The mean and standard deviation for the samples including the upper 
limits are shown with the red squares and error bars (see Table [5J . We also compare the mean and standard deviation of the complete 
early-type galaxies (S0+S0a and E) and the complete sample of spiral galaxies (C12) in blue. The stellar masses are from L. Cortese et al. 
(in preparation). 
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Figure 10. Top panel: the FIR luminosity from the cold dust component Lpipt(cold) vs. the B-band luminosity (left). Blue points arc 
SO+SOa galaxies and red points are the ellipticals. Upper limits are plotted as downward-pointing triangles. The gray solid line shows the 
expected relationship (arbitrarily normalized) for Lpia oc Lg 1 if the FIR emission is produced by dust from stellar atmospheres. The best-fit 
line to the SO+SOa sample including the upper limits is shown by the dashed blue line (m = 1.22, c = —3.53, o> = 0.39) as estimated from 
the Buckley James method. The correlation parameters for ellipticals are shown by the dashed red line (m = 1.00, c = —2.44, oy = 0.3). 
Right: the FIR emission vs. X-ray luminosity {Lx) normalized by the blue luminosity (Lb). Bottom panel: the dust mass M d vs. the B- 
band luminosity (left) and vs. X-ray luminosity (right). The logM^ vs logLg relationship is described with m = 1.11, c = — 5.08,o> = 0.58 
for the SO+SOa (blue dashed) and m = 1.89, c = —14.40, oy = 0.97 for the ellipticals (red dashed). The dashed gray lines represent the 
loci where dust is replenished by stellar mass loss and destroyed by dust sputtering in the hot gas (the labels represent the log of the dust 
destruction timescale assumed). This model assumes a gas-to-dust ratio of 150. The errors on the dust mass are similar to the symbol size 
for most galaxies, and a representative error bar is shown. 
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Table 6 

Results of the Two-sample Tests 



Sample 


Parameter 


P(GGW) 


P(LR) 


P(PPGW) 


SO X E 


M d 


0.0005 


0.0006 


0.004 




M d /M t 


0.000 


0.000 


0.000 




M(H 2 ) 


0.007 


0.009 


0.003 


Virgo X not Virgo 


M d 


0.14 


0.09 


0.14 




M* 


0.189 


0.151 


0.165 




M d /M« 


0.055 


0.082 


0.078 




M(Hi) 


0.008 


0.011 


0.005 



Notes. — The S0+S0a versus elliptical (top) and Virgo versus non- 
Virgo galaxy samples tested for the probability they are drawn from 
the same population. The tests include GGW: Gehan's Generalized 
Wilcoxon test; LR: log-rank test; PPGW: Peto and Peto Generalized 
Wilcoxon test. These are determined using the IRAF STSDAS task 

TWOSAMPT. 

2010; Serra et al. 2009; di Serego Alighieri et al. 2007) 
and detections of molecular gas have been made in many 
(Sage et al. 2007; Lucero <fc Young 2007; Combes et 
al. 2007). A number of studies indicate ongoing star 
formation in a significant fraction of ETGs (Kaneda et 
al. 2005, 2008; Bressan et al. 2006; Panuzzo et al. 2007, 
2011; Temi et al. 2007a, 2007b) and the detection of dust 
in 50% of our sample strongly supports the conclusion 
that a significant fraction of ETGs contain a cool ISM. 

We can compare the ISM in the ETGs with that in late 
types in two other ways. First, using the measurements 
of the dust, atomic, and molecular gas that exist for eight 
galaxies from our sample (Table |4|), we estimate that 
the mean gas-to-dust ratio is io 2 08±0 08 ; similar to the 
typical value for late-type galaxies: ~ 100-200. 

Second, in late-type galaxies there is a tight correlation 
between FIR and radio emission (Wrobel & Heeschen 
1988; Helou et al. 1985; Devereux & Eales 1989). This is 
usually explained as the result of a correlation between 
the FIR emission and the SFR, with FIR emission origi- 
nating either from dust heated by young stars (e.g., Vla- 
hakis et al. 2007 and references therein; G. P. Ford et al. 
in preparation) or due to dust itself tracing the gas mass, 
which in turn fuels star formation (Rowan-Robinson et 
al. 2010; Bendo et al. 2011). The radio emission is then 
also correlated with the SFR as it originates from rela- 
tivistic electrons produced by supernova remnants after 
the young stars die. Whether or not this correlation is 
truly caused by the SFR, we can at least see whether 
the ETGs in our sample have a similar ratio of FIR- 
radio as spiral galaxies. Sixteen ETGs are detected at 
1.4 GHz (see Figures [T] and [2]) by the Faint Images of the 
Radio Sky at Twenty centimeters (FIRST) radio survey 
(Becker et al. 1995). For these sources we have estimated 
the value of the parameter, q, introduced by Helou et al. 
(1986), defined as 



log 



L 



FIR 



3.75 x 10 12 Wnr 



log 



Li 



4 GHz 



Wm 



-2 



where Lfir and ^i.4ghz are the FIR and radio lumi- 
nosities, respectively. Yun et al. (2001) found a median 
value for nearby galaxies of q = 2.64 ± 0.02, with most 
star- forming galaxies having values of q between 2 and 3. 
Wrobel & Heeschen (1988), and more recently Combes 
et al. (2007) and Lucero & Young (2007), have found 
similar values of q for ETGs in which star formation is 



occurring. 

Two of our galaxies (NGC4636, NGC4374) have q < 
1.8, which suggests they host radio-loud AGNs. Thir- 
teen out of 16 sources with FIRST detections in our HRS 
sample have 2.15 < q < 3.32, similar to star-forming late 
types. Nine of these sources also have CO detections 
(Table 0]) and therefore 70% of the galaxies which the q 
values suggestive of star formation also have a reservoir 
of molecular gas. These galaxies lie below the 'quiescent' 
UV-optical boundary defined as NUV-r < 5.4 (Section 
13.21) and a literature search reveals evidence of residual 
star formation in seven out of the nine galaxies, with sig- 
natures including mid-IR emission (Panuzzo et al. 2007; 
Shapiro et al. 2010), line diagnostics (Sil'chenko et al. 
2010; Crocker et al. 2011), and UV emission (Cortese & 
Hughes 2009). 

Converting the q ratio into SFRs for ETGs is extremely 
complex due to the difficulty in disentangling the AGN 
component, other processes that may be heating the 
dust, and the thermal contribution to the radio emis- 
sion. If the q ~ 2 ratios for those galaxies with molec- 
ular gas are an indication of star formation, then there 
are not only cool gas and dust in these galaxies but also 
stars forming at a measurable rate for at least 15% of 
the sample. We defer a full analysis of SFRs for future 
work, though we note that a literature search reveals 
SFRs ranging from 0.03 to 0.40 Mgyr" 1 in 23% of the 
HRS ETG sample (see references listed above). 

In summary, although the mass of the ISM in ETGs is 
less than the mass of the ISM in late types, observations 
of the radio continuum, CO, and Hi are consistent with 
the ISM in ETGs being quite similar to that found in 
late types. The only difference is that, on average, the 
dust in ETGs has a temperature that is a few degrees 
higher, possibly due to the dust grains in early types 
being exposed to a more intense interstellar radiation 
field or from the more centrally distributed dust emission 
in ETGs (e.g., Sauvage & Thuan 1992; G. Bendo et al. 
2011, submitted). 

5. DISCUSSION 

5.1. Evidence for dust-depleted disks in SOs 

The large fall in the dust-to-stellar-mass ratio between 
early- type spirals and S0+S0a galaxies (Figures [8] and [9]) 
suggests that SOs contain a smaller mass of ISM per mass 
of stars than early- type spirals (i.e., those with Hubble 
type ranging from Sa to Sbc). However, the bulge-to- 
total mass ratio is known to be larger for SOs than early- 
type spirals: is the fall in the dust-to-stellar-mass ratio 
we observe here simply a reflection of the more massive 
bulges of SOs? We can test this idea by making the as- 
sumptions that (1) bulges do not contain dust and (2) the 
dust-to-stellar-mass ratio is the same for all disks. With 
these assumptions, to explain the factor of 10 decrease in 
the dust-to-stellar-mass ratio of SOs relative to early-type 
spirals, the bulge-to-total stellar-mass ratio must change 
by a very large factor. If the bulge-to-total stellar-mass 
ratio is 0.1 for early-type spirals, it must be 0.91 for SOs. 
For larger values of the bulge-to-disk stellar-mass ratio 
for the early- type spirals (Sa-Sc), even larger ones are 
required for the SOs. Weinzirl et al. (2009) have found 
that the proportion of the total mass of a galaxy that 
is in the bulge increases from ~0.1 for Sbc galaxies to 
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Figure 11. Ratio of the submm and optical sizes for the HRS 
sample. The red histogram shows the distribution of Dsubmm/^opt 
for the early-type spiral galaxies (Sa-Sbc) with the same color cut 
as the SO and ellipticals, i.e., K < 8.7 (the original sample of early 
types is presented in Cortese et al. 2010b and C12). The green 
shaded histogram shows the distribution for the SOs in this paper, 
with the blue shaded histogram showing the SOs that are not in one 
of the Virgo clouds or the outskirts of Virgo (Boselli et al. 2010b). 

^0.2 for SO+SOa galaxies (their Figure 14). In larger 
samples, Graham & Worley (2008) and Laurikainen et 
al. (2010) find values for this ratio of approximately 0.1- 
0.2 for Sab-Sbc galaxies and 0.3-0.5 for SO galaxies (their 
Figure 4), a somewhat larger change but still too low to 
explain the different dust-to-stellar-mass ratios. 

Another way to investigate this is to compare the ra- 
tio of the submm diameter to the optical diameter (e.g., 
Cortese et al. 2010b), where the optical diameter is likely 
to be a good measure of the size of the disk. We mea- 
sured this ratio for two samples: (1) all of the SO galax- 
ies and (2) all early-type spiral galaxies in the HRS with 
K < 8.7 (the limit used for ETGs). For the optical 
size of the galaxy we used ma j O r, the diameter along 
the major axis measured to the standard B 25th mag 
arcsec -2 isophote (£>(25)). For the submm size, we used 
the same definition, finding the ellipse that provided the 
best fit on the 250 /im image to an isophotal brightness of 
6.7 x 10~ 5 Jy arcsec~ 2 . Figure [TT] shows the distributions 
of D su b mm /D opt for the S0+S0a galaxies and the early- 
type spirals with the same .ftT-band magnitude selection; 
these populations are significantly different (P < 0.001 
leveP 7 !) . We do not apply a correction for beam smear- 
ing, which increases the angular sizes of the objects, an 
effect which is biggest for the objects with smallest an- 
gular size. If a correction was made for this bias, the 
difference between SO+SOa and early-type spiral popu- 
lations seen in Figure[TT] would increase. We have also 
compared the submm-optical sizes for the SO galaxies 
outside of Virgo in Figure QTJ The values for these are 
still clearly different from the early-type spirals, strongly 
suggesting that the dust-depleted disks of SOs are not 
caused by a current cluster-environmental effect. 

27 Using the Kolmogorov-Smirnov two-sample test. 



Both approaches imply that the dust-to-stellar-mass 
ratio is lower for the disks of SO galaxies than the disks 
of early-type spirals, suggesting that the disk of an SO 
contains a smaller mass of ISM than the disk of an early- 
type spiral with the same stellar mass. This is not an 
entirely surprising result, since for the last 40 years the 
working definition of an SO galaxy is that it is a galaxy 
with a disk but no sign of spiral arms or, if it is an edge-on 
galaxy, of a dust lane (Sandage 1961). Nevertheless, even 
if the qualitative result is not surprising, Herschel has 
allowed us to determine how little ISM the disks of SO 
galaxies do contain. 

5.2. Ellipticals and SOs or slow rotators and fast 
rotators ? 

So far in this paper, we have classified ETGs based on 
their optical morphology into ellipticals and SOs. Based 
on the ATLAS 3D study of the stellar kinematics of 260 
ETGs, Cappellari et al. (2011) have argued that a 
physically more meaningful way is to divide ETGs into 
slow rotators and fast rotators. Emsellem et al. (2011) 
found that 66% of the galaxies in ATLAS 3D , which 
are traditionally morphologically classified as ellipticals, 
have disk-like kinematics and are therefore fast rotators. 
Given the similarity in the kinematics, the suggestion 
is then that early-type fast rotators are part of the same 
evolutionary sequence as the late types, but arc at differ- 
ent stages of their evolution. Slow rotators are defined as 
galaxies with bulge-like properties, likely to be the end- 
point of systems that have undergone complex merger 
histories; these are massive, lie on the red sequence, and 
contain very little cool interstellar material. 

Eleven (18%) of the ETGs in our sample are slow rota- 
tors, including nine ellipticals (39% of the ellipticals) and 
two SOs (5%). In Figure [TW a). we plot dust-to-stellar- 
mass ratio versus the quantity Afl c/2 /^/e e /2, where e is 
the ellipticity of the galaxy and R e is the effective ra- 
dius. This parameter is suggested by the ATLAS 3D team 
as appropriate to separate out the slow and fast rotators 
since it is a good measure of the projected stellar angular 
momentum per unit mass (Emsellem et al. 2007). There 
is no sign of a correlation between the two parameters for 
galaxies in our sample. Indeed, with Herschel we have 
detected dust emission from 50% of the slow rotators, in- 
cluding four ellipticals, which is not very different from 
the detection fraction for all the ETGs (Table . 

The fact that there is a clearer distinction in the dust 
properties of the ETGs when they are put in morpho- 
logical classes (Figures [H] and ^ than when they are put 
into kinematic classes does not mean that the kinematic 
division is not physically meaningful. It is possible that 
the morphological classification is more connected to the 
mass of ISM in a galaxy than the kinematic classifica- 
tion, although the kinematic definition may be telling us 
more about the evolutionary state of the galaxy. As Her- 
schel detects a number of slow rotators and the lack of 
correlation seen in Figure I12f a) , this may be additional 
evidence that the cool ISM in ETGs is acquired by ran- 
dom gravitational encounters. Again, the amount of dust 
may not be connected with the evolutionary state of the 
galaxy. In Figure I12f b) , we show the same diagram but 
this time labeling the galaxies by their Hi mass. There is 
a tendency for galaxies with larger atomic gas masses to 
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Figure 12. Top: the dust-to-stellar-mass ratio vs. the kinemati- 



cal parameter ^R e /2/-. 



where e e /2 is the ellipticity measured 



within an aperture defined at 0.5H e and R e is the effective ra- 
dius of the galaxy. \r is a proxy for the angular momentum per 
unit stellar mass and the ^/e e /2 term accounts for the shape of a 
galaxy since a more flattened galaxy would be expected to have 
a stronger anisotropy. Blue points are S0+S0a galaxies and red 
points are ellipticals. The gray solid line marks the region defined 
as the boundary between fast and slow rotators (Emsellem et al. 
2011), i.e., between galaxies with disk-like rotation and bulge-like 
kinematics. Bottom: as above but color-coded with respect to the 
Hi mass. Galaxies which are known to be mergers or have obvious 
dust lanes in the optical are indicated with + (mergers) and X (dust 
lanes). Galaxies without filled circles are the Hersc/ieZ-detected 
galaxies that have no corresponding Hi observations. Many of the 
galaxies at the lower end of the Hi mass scale are upper limits 
(Table [4t . The gray dashed line roughly marks the detection limit 
of the survey and note that NGC4486 (HRS183) is not shown in 
these plots. 

have higher dust-to-stellar mass ratios. Since the 21 cm 
morphologies have often been adduced as evidence for an 
external origin of the ISM (e.g., Serra et al. 2009), this 
figure may be additional evidence that the mass of the 
ISM is rather unconnected to the evolutionary history of 
the stellar populations in the galaxies. 

Finally, we note that connections between the dust 
properties and stellar kinematics of the ETGs may be- 
come evident in a larger sample. The HRS sample of 
ETGs is a statistically complete quasi-volume-limited 
sample, but the volume contains the rich environment 
of the Virgo Cluster and is therefore not a representative 
sample of the local universe. The HRS ETGs also cover 



a limited range of stellar mass (Section |31 Figure [3]). An 
important future project would be to observe the dust 
in a sample of ETGs containing more galaxies in low- 
density environments and with a wider range of stellar 
mass. 

5.3. Galaxy evolution 

We can use the results from the Herschel observations 
to make an inference about the evolution of early types. 
On the assumption that the dust we observe is delivered 
to the galaxies externally, we can make a rough estimate 
of the effect such an interaction would have on the stel- 
lar mass of the ETG. We use the mean values of M^/M* 
shown in Figure [5] (Table [5]) for the different morpho- 
logical classes, (10~ 5,9 for ellipticals, 10~ 44 for SOs, and 
~ 10~ 30 for spirals). On the assumption that all the 
dust we detect in an ETG is the result of a past merger 
with a spiral, the ratio of the mass of the spiral to the 
original ETG must be ~ 0.05 for an elliptical and ^0.1 
for an SO. These estimates suggest that such interactions 
are only minor mergers and do not represent significant 
galaxy-building events. The estimates are also upper lim- 
its because the dust may have been acquired by a tidal 
interaction without the galaxies necessarily merging, and 
some of the dust may also have been produced inside the 
galaxy. 

6. CONCLUSIONS 

We present the FIR and submm observations with 
Herschel of 62 ETGs, including 23 ellipticals and 39 
S0+S0a galaxies. We find the following results: 

1. We detect 24% of the ellipticals and 62% of the 
SOs. The optical and X-ray luminosities of the 
detected sources show we are detecting the most 
massive ellipticals, though the SO sample is rep- 
resentative of the general population. Of the ten 
'pure' ellipticals in our sample only one is detected 
by Herschel (M86), while six of the eleven 'impure' 
ellipticals are detected. The detection rate for the 
ETGs outside Virgo appears to be higher than for 
those inside the cluster but this is not statistically 
significant. 

2. The mean dust masses for the detected galaxies 
are logM d = 6.3 ± 0.1 and 5.5 ± 0.3 M for the 
SO+SOa and E population, with average dust-to- 
stellar-mass ratios of log(M t ;/M st ) = — 4.1±0.1 and 
—5.3 ± 0.2. Including the upper limits, the aver- 
age dust masses for the detected and undetected 
SO+SOa and elliptical sources are logAf t j(M ) = 
5.9 ±0.1 and 5.2 ±0.1, with dust-to-stellar-mass 
ratios of log(M d /M*) = -4.4 ±0.1 and -5.8 ±0.1. 

3. The mean dust temperature for the ETGs de- 
tected by Herschel is ~ 24 K, warmer than the dust 
in late-type galaxies. The entire early-type sam- 
ple, including non-detections, has a dust mass of 
logM d = 5.6 ± 0.1 and log(M d /M*) = -5.1 ± 0.1. 
In comparison, the average dust-to-stellar mass ra- 
tio for spiral galaxies in the HRS is —2.59 ± 0.03. 
The dispersion in the dust-to-stellar mass ratio is 
much greater for ETGs than spirals. 
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4. The NUV-r colors show that virtually all the ETGs 
lie close to the red sequence, but there is evidence 
from UV and optical colors, radio continuum ob- 
servations, and literature searches that a ~ 20% of 
the sample have had a recent star-forming epoch or 
have significant residual star formation. However, 
the NUV colors have not been corrected for extinc- 
tion or contamination from the old stellar popula- 
tion. Our ETGs are redder in NUV-r with lower 
dust-to-stellar mass ratios than the ETGs detected 
in the H- ATLAS survey. 

5. We show that the detection of cold dust, the ratio 
of far-infrared to radio emission, and the gas-to- 
dust ratios all indicate that many ETGs contain a 
cool ISM in which stars are forming, similar to that 
seen in spiral galaxies. 

6. We find no evidence for a correlation between the 
FIR luminosity or dust mass of ETGs with the op- 
tical luminosity, suggesting that the main source of 
the dust in the galaxies detected by Herschel is not 
mass loss from evolved stars. This, together with 
the large spread in the dust-to-stellar-mass ratio of 
ETGs, suggests a significant fraction of the dust in 
these galaxies is acquired externally via mergers or 
tidal interactions. 

7. We use the results from the ATLAS 3D survey to 
divide our sample into fast and slow rotators. We 
show that the difference in dust properties between 
SOs and ellipticals is more obvious than the differ- 
ence between slow and fast rotators. We suggest 
that this may be additional evidence that the dust 
in ETGs has been acquired by gravitational en- 
counters and may not tell us much about the evo- 
lutionary state of the galaxy. 

8. The low dust-to-stellar-mass ratios of SOs com- 
pared with early-type spirals cannot be explained 
by the larger bulge-to-disk ratios of SOs. The rela- 
tive sizes of the dust sources in SOs are also smaller 
than seen in early- type spirals. These suggest that 
the disks in SOs contain much less dust (and pre- 
sumably gas) than spiral disks with a similar size. 
This effect is probably not being caused by cur- 
rent environmental processes adding to the weight 
of evidence that early-type spirals are not being 
transformed into SOs in significant numbers at the 
current epoch. 

9. If the cool ISM in ETGs is acquired as the result 
of a tidal interaction or merger, an upper limit on 
the increase in the stellar mass of the ETG due to 
the interaction is 1% for the ellipticals and 10% for 
the SOs, suggesting that the interactions are not 
significant galaxy-building events. 

Our sample of ETGs has lower dust-to-stellar-mass ra- 
tios and dust masses compared to previously published 
Herschel samples. The results from our sample are im- 
portant for interpreting chemical evolution models which 
model the production and destruction of dust in ellipti- 
cals as well as providing a low-redshift benchmark for un- 
derstanding the evolution of dust and gas in high redshift 



surveys (e.g., Rowlands et al. 2012, Dunne et al. 2011). 
Finally, we note that to obtain a better understanding 
of some of these issues it will be important to follow-up 
our statistical study with detailed observational studies 
of individual ETGs, for example with ALMA. 

We thank everyone involved with the Herschel Ob- 
servatory. PACS has been developed by a consor- 
tium of institutes led by MPE (Germany) and in- 
cluding UVIE (Austria); KU Leuven, CSL, IMEC 
(Belgium); CEA, LAM (France); MPIA (Germany); 
INAF- IFSI/OAA/OAP/OAT, LENS, SISSA (Italy); 
and IAC (Spain). This development has been sup- 
ported by the funding agencies BMVIT (Austria), ESA- 
PRODEX (Belgium), CEA/CNES (France), DLR (Ger- 
many), ASI/INAF (Italy), and CICYT/MCYT (Spain). 
SPIRE has been developed by a consortium of institutes 
led by Cardiff University (UK) and including University 
of Lethbridge (Canada); NAOC (China); CEA, OAMP 
(France); IFSI, University of Padua (Italy); IAC (Spain); 
Stockholm Observatory (Sweden); Imperial College Lon- 
don, RAL, UCL-MSSL, UKATC, University of Sussex 
(UK); and Caltech/JPL, IPAC, University of Colorado 
(USA). This development has been supported by na- 
tional funding agencies: CSA (Canada); NAOC (China); 
CEA, CNES, CNRS (France); ASI (Italy); MCINN 
(Spain); Stockholm Observatory (Sweden); STFC and 
UKSA (UK); and NASA (USA). HIPE is a joint devel- 
opment by the Herschel Science Ground Segment Con- 
sortium, consisting of ESA, the NASA Herschel Sci- 
ence Center and the HIFI, PACS and SPIRE consor- 
tia. This research made use of the NASA Extragalac- 
tic Database (NED: |http://ned4 pac.caltech.edu/| which 
is operated by the Jet Propulsion Laboratory, Califor- 
nia Institute of Technology, under contract with the 
National Aeronautics and Space Administration. The 
research leading to these results has received fund- 
ing from the European Community's Seventh Frame- 
work Programme (/FP7/2007-2013/) under grant agree- 
ment no. 229517. This research made use of APLpy, 
an open-source plotting package for Python hosted at 
http://aplpy.github.com. We thank Edward Gomez for 
useful and informative discussions and H.L.G. acknowl- 
edges the support of Las Cumbres Observatory. 

Facilities: Herschel (PACS and SPIRE), Spitzer 
(MIPS), IRAS 

APPENDIX 

A. NOTES ON THE ELLIPTICALS DETECTED BY 
HERSCHEL 

HRS3, NGCS226. A semicircular lane of dust in 
NGC3226 is seen as absorption in the optical (Martel 
et al. 2004) with dusty strands extending north-south. 
The source is known to be interacting with the nearby 
spiral NGC3227 (the bright source to the south of the 
elliptical — Figure [2J, which has clear tidal gas trails 
(Mundell et al. 2004). The elliptical appears slightly 
extended in the 250 /jm image but more compact than 
in the optical image (Figure [2]) • NGC3226 is part of the 
Leo Cloud. 

HRS138, NGC4374 (M84). M84 has a radio bright core 
and radio jets. Hubble Space Telescope (HST) images 
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show dust lanes across the center of the galaxy (Bower 
et al. 1997), which are approximately perpendicular 
to the radio jet. The galaxy was detected at 850 fim 
with SCUBA by Leeuw et al. (2000). They concluded 
that most of the 850 /im flux was synchrotron emission, 
but that the emission at shorter wavelengths detected 
by IRAS was produced by ~ 10 5 M© of dust at a 
temperature of 35 K (see also Boselli et al 2010a). The 
250 /im source is unresolved and is coincident with the 
radio core (Figur^2]). NGC4374 is a member of the 
Virgo Cluster. 

HRS150, NGC4406 (M86). M86 is a well-known IR- 
bright giant elliptical. Two dust features were detected 
with IRAS and originally attributed to dust stripped 
from M86 due to its motion through the cluster. The 
discovery of atomic gas offset from the center of M86 
and decoupled from its stellar disk supports a tidal 
interaction (Li & van Gorkom 2001), confirmed when 
Kenney et al. (2008) detected strong Ha features 
extending from M86 to the nearby spiral NGC4438. The 
distribution and velocity of the ionized gas provide clear 
evidence for tidal interaction between these two giants. 
Herschel SPIRE observations of M86 and NGC4438 
showed that the dust emission is spatially correlated 
with the ionized gas between the two galaxies (Gomez 
et al. 2010; Cortese et al. 2010a) implying that the dust 
is material stripped from the nearby spiral. In Figure [21 
the appearance of the dust emission in M86 is clearly 
different to the other ellipticals in our sample, with faint 
filamentary features seen within D(25). NGC4406 is a 
member of the Virgo Cluster. 

HRS183, NGC4486 (M87). M87 is the brightest galaxy 
in the Virgo Cluster and well known for the jet extend- 
ing from the nucleus seen at radio, optical, and X-ray 
wavelengths (Junor et al. 1999). The presence of dust 
was inferred from HST g - z color images (Ferrarese 
et al. 2006) and the strong FIR emission seen with 
Spitzer (Perlman et al. 2007). The latter result might 
be explained by synchrotron emission from the central 
radio source. Although the optical images show dust is 
present, Herschel observations of M87 with PACS and 
SPIRE as part of HeViCS (Baes et al. 2010; Boselli 
et al. 2010a) found no evidence of a dust component 
in excess of the synchrotron emission in the FIR and 
submm. Baes et al. place an upper limit on the dust 
mass of 10 5 Mq. In the Herschel image, M87 is a bright, 
extended source. 

HRS186, NGC4494. NGC4494 is often described as an 
'ordinary elliptical' (Capaccioli et al. 1992). O'Sullivan 
& Ponman (2004) found that its X-ray luminosity was 
two order of magnitudes lower than expected for its 
optical luminosity and there are signs of a small dust 
disk in absorption (Tran et al. 2001). X-ray faint galax- 
ies such as NGC4494 may arise due to losing their hot 
X-ray gas in outflows. Further evidence of interactions 
with nearby galaxies or mergers is the low metallicity 
estimated from the X-ray gas (< 0.1Z Q , O'Sullivan 
& Ponman 2004), and this may indicate dilution of 
interstellar material via an infall of unenriched (cold) 
material. The galaxy is slightly extended at 250 /im. 
NGC4494 is a member of the Coma I Cloud. 



HRS241, NGC4636. NGC4636 has dust features and 
an unusual X-ray morphology (Temi et al. 2007b). 
The origin of this morphology is thought to be recent 
outbursts from the central AGN (Jones et al. 2002). 
Temi et al. (2003) used ISO observations to show that 
the dust mass for NGC4636 is far in excess of that 
expected from stellar mass loss and proposed that the 
dust was accreted in a very recent merger with a dusty, 
gas-rich galaxy (similar to M86, Gomez et al. 2010). 
NGC4636 is part of the Virgo Cluster. 

HRS258, NGC4697. NGC4697 is an X-ray-faint galaxy 
given its optical luminosity, which may be a result of 
severe loss of interstellar gas via stripping or outflows 
(Sarazin et al. 2001). As in NGC4494, the metallicity of 
the gas is low (< 0.07Z Q ), pointing toward dilution of 
the interstellar medium with unenriched cold gas. The 
source is slightly extended at 250 /im. NGC4697 is in 
the outer regions of the Virgo Cluster. 

B. FIR EMISSION VERSUS OPTICAL AND X-RAY 

We investigated whether there is a correlation between 
Lx and Lb using the statistics in the ASURV package 
for dealing with censored data. We find evidence for a 
correlation for the elliptical galaxies (P = 99.390. The 
relationship between the two parameters for ellipticalJ^l 
is described by the following relationship with deviation 
from the regression cr r = 0.43: 

\ogL B = (2.95 ± 1.03)logi x + 10.79 (Bl) 

This relationship agrees with other IR samples domi- 
nated by massive ETGs (see also Brown & Bregman 
1998), whereas X-ray studies indicative of the entire pop- 
ulation find are described by Lx oc L B 3 (O'Sullivan et 
al. 2001; Temi et al. 2004). 

If the dust produced in ETGs is due to stellar mass 
loss, the dust mass would be roughly proportional to the 
mass of the stars, and hence to Lb- If the dust is then 
widely dispersed and heated by the starlight as well as 
electron collisions in the hot, X-ray emitting gas, the FIR 
emission will also depend on the density of the gas or on 
the stellar density, which are both roughly proportional 

to L 1 ^ 2 . Therefore, we would expect Lfir oc L B A ~ l e . 
Bregman et al. (1998) indeed found this relationship was 
true for JiMS-detected ETGs, though Temi et al. 2004 
found no correlation for their /S'O-detected sources. In 
the Herschel sample, we see that the FIR luminosity of 
SO+SOa galaxies does appear to increase with the optical 
luminosity, although this is not a statistical correlation, 
with Spearman ranlFl coefficient r$ = 0.36 and prob- 

28 Using the Kendall r test, a statistic used to measure the 
association between X and Y; this test is appropriate for small 
samples with upper limits where the underlying distributions of X 
and Y are not known. 

29 Using the Buckley- James method, a standard linear regression 
estimator. This test requires that the censoring distribution about 
the fitted line is random and does not require the residuals to be 
Gaussian as other regression methods do. 

30 Using the Spearman rho correlation test; this statistic is more 
appropriate than the Kendall r used for the ellipticals since the 
SOs are a larger sample with more detections. The correlation 
determines how well the relationship between X and Y can be 
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ability that a correlation is present of P = 97%. The 
SO+SOa data can be fit by the following equation^ with 
cr r = 0.39: 

log£ F rR,(S0) = (1.22 ± 0.42)log (L B ) - 3.53. (B2) 

The FIR luminosity of the elliptical galaxies increases 
with Lb but again, we find that the correlation is not sta- 
tistically significant, with Lb (t = 0.36 and P = 91%P), 
the results of the regression analysis produces the follow- 
ing relationship (with tr r = 0.3323): 

logi F m(E) = (1.00 ± 0.77)log (L B ) - 2.44. (B3) 

The correlations in Equations (|B2|) and (|B3|) are plot- 
ted in Figure [TUJ The lack of any strong correlation be- 
tween Lfir and L b does not provide any evidence for the 
hypothesis that the dust responsible for the FIR-submm 
emission is produced by stellar mass loss, although the 
ETGs do at least fall in roughly the right place in the 
figure. 

In Figure [TU] (lower panel), we also compare the dust 
mass with optical and X-ray luminosities but find no ev- 
idence for a correlation due to the small numbers in the 
samples with both Herschel detections and X-ray fluxes 
in the literature. 

There is also no evidence for a significant correlation 
between dust mass and Lb for SO galaxies (rg = 0.27 
and P = 90%), with a best-fit relationship (a r = 0.5&P 3 ): 

logM rf (S0) = (1.11 ± 0.62)log£ B - 5.08 (B4) 

The dust mass and optical luminosity for the ellipticals 
are also not correlated (r = 0.30 and P = 90%^), with 
a best-fit relationship (cr r = 0.87^ 

logAfd(E) = (1.89 ± 1.99)logL B - 14.40 (B5) 

The remaining panel in Figur dTUl shows dust mass 
rather than FIR luminosity plotted against X-ray lumi- 
nosity, with both quantities being normalized by Lb! the 
correlation is weaker when FIR luminosity is converted 
into dust mass. 
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